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ABSTRACT 


The present work deals with the steady state analysis and 
control of torque harmonics of the induction motor fed by three 
phase current source inverter. 

The analytical solutions for the rotor currents, stator 
voltages and torque waveform in the stationary dq frame 
attached to thestator have been obtained. The rotor currents 
and torque harmonic spectrums have also been obtained through 
frequency domain analysis. 

The next part of the thesis deals with tho control of 
the torque harmonics by modulation of dc input current. A 
detailed study has been made of the modulation by exponential 
and cosinusoidal v/aveforms and their effects on various torque 
harmonics. A general nonlinear second order equation has been 
derived for obtaining the waveshape of the inverter input 
current required to produce a desired torque waveform. A 
general analytical solution of this equation has been done for 
the case of stationary rotor. Further for the case of rotating 
rotor and constant torque, this equation has been solved using 
phase plane analysis and numerical integration. The current 
profile for producing constant torque at any rotor speed in tho 
motoring mode has been obtained. 



CHAPTER 1 


INTRODUCTION 

It is well known that an induction motor supplied by a 
current source inverter draws line currents with an approxi- 
mately quasi square waveform [l]. Consequently, the general 
analysis of the motor with sinusoidal voltages and currents is 
not directly applicable in this case. The most popular 
approaches for the analysis of ac motor drives, which employ 
non sinusoidal sources are harmonic superposition techniques 
using either symmetrical components [2] or method of multiple 
reference frames [s]. The analysis of the performance of the 
induction machine fed by a current source inverter has been 
carried out in [l] , [4] , [5] , [6] , [?] using dq axes model of 
induction motor. Several methods of solving dq equations 
have been proposed. In [4] , [5] , [6] , the analysis of induction 
motor is carried out assuming ideal waveshape for stator current 
waveform. In [7] the analysis has been done using digital simu- 
lation. None of these give direct analytical expressions for 
calculating referred rotor currents stator voltages and torque 
waveform. In this thesis analytical expressions are obtained 
using d-q model of the machine for referred rotor currents, 
stator voltages and electromagnetic torque. Analysis has been 
done for both the cases of ideal current source inverter and 
the inverter with nonzero commutation time. For the latter, 
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therisc-' and fall of the stator currant during 
been taken of cosinusoidal nature [l]. 


commutation has 


It is possible to control the torque harmonics of an 
induction motor fed by a current source inverter either by 
modulating the dc current input of the inverter or by modula-" 
tion within the inverter. In [s] , [9] the effects on the 
torque harmonics by modulating the inverter input current by a 
particular exponential v^aveform has been studied. Modulation 
ivithin the inverter is referred to as pulse vddth modulation in 
literature [lO], Reference [lO] gives the effects on torque 
harmonics because of tho pulse width modulated current source 
invoDrters. In this thesis a detailed study has been made of 
che mociulation of inverter input current with exponontlal and 
cosinusoidal waveforms. It has been shown that the parameters , 
of the modulating waveforms control the torque harmonics. 
Optimum parameters for eliminating the dominant sixth harmonic 
have been calculated. A general nonlinear second order equa“ 
xron has also been derived for obtaining tho v^avoshapo of the 
invertor input current required to produce a desired torque 
waver orii. A general analytical solution of this equation has 
been done for the case of stationary rotor. Further for tho 
case of rotating rotor ’and constant torque this equation has 
been solved using phase piano analysis and numerical integra- 


tion. 



jThe c’q model of the induction motor, with axes attached 
to the str>tor has been used throughout for the analysis. A 
brief description of dq axes model for three phase induction 
■aotor is given in Sec. 1,1. 

In Chagtc-r_2, the steady state analysis of three phase 
induction motor fed by an ideal th:ceo phase current source 
inverter has been done to compute the referred rotor current 
and the stator voltage waveforms. Tho analytical solutions are 
detained for rotor currents and stator voltage for three inter-^ 
vals of the inverter period. Boundary conditions are verified 
analytically for the c'’se of stationary rotor. The resulting 
vxoressions for the case of the rotating rotor are complicated 
and hence, closed form equations have not been obtained. A 
conputer program has boon wD'itten to compute tho rotor currents, 
/.naiyticnl solution results in a considerable reduction of com- 
puter time compared to the numerical solution of dq equations. 
Tho i'<oundary conditions for the case of rotating rotor arc 
verified with the computer program. The values of the rotor 
current harmonics have also been obtained through frequency 
domain analysis of the induction motor. It has been shown that 
the rotor harmonics as obtained from tho analytical solution 
ore in close agreement to those obtained through frequency 
domain analysis. 

In Chapter 3, tho analysis similar to that in Chapter 2 
has been done. In this case tho invertor is assumed to have 
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a nonzero commutation time. During this commutation time^ 
too rise and fall of the current is approximated by a 
consinusoidal function [l]. 

In Chapter 4__o the electromagnetic torque of the current 

fed induction motr;r has been computed. Two methods to compute 

tl'a'i torque have been given. One obtains the analytical solu- 

ti;.'*'' of torque using the time domain expressions of stator and 

referred rotor currents in dq frame. The other computes the 

torquo harmonic spectrum using the harmonic components of the 

stator and rotor currents in d-q frame. Using these procedures- 

the torque- waveform and spectrum has been studied for the case 

of ideal current source invertor supplying an induction motor, 
to the 'adulation of the inverter input dc cc’crent by '-^xponon"- 

tial r;e'-iulv^ tion. 

In Cha ’ter__5, a detailed study has been made of the 
effects on torquo harmonics duo to the modulation of the inver- 
tor input dc current by exponential and cosinusoidal waveforms. 
The exponential modulation for stationary rotor case gives con- 
stant torque. Cosinusoidal type of waveform may be present in 
the invertor input current due to imperfect filtering in the 
dc link. It has been shown that it is not possible to get an 
exactly constant torque with these* typos of modulation for'jthe 
case of rotating rotor. However, by a certain choice of the 
paramotors of the modulating waveforms, it is possible to 
reduce certain torque harmonic. 
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Chapter 6 deals with the calculation of stator current 
profile of the inrluction motor fed by current source inverter 
to produce a specified torque waveform. A general nonlinear 
second order equation has been derived for obtaining the 
v/avoshapo of the inverter input current required to piroduce 
ti'io desired torque vravoform. A general analytical solution of 
this equation has boon done for the caso of stationary rotor. 
Further for the caso of rotating rotor and constant torque 
this equation has boon solved using phase plane analysis and 
numerical integration. The current profile for obtaining a 
constant torque at any rotor speed in the motoring mode has 
been obtained. 

1.1 DQ MODEL FOR INDUCTION MOTOR 

The development of the dq model for the induction motor 
using arbitrary revolving frame has boon done by Krause [s]. 
In this section the equations of motor performance and the 
transformations from throe phase balanced system to the dq 
frame have boon outlined briefly using [3]. 

Fig. 1.1 shows the angular relation of the stator and 
rotor axis of a throe phase machine with the third set which 
is an orthogonal sot (dq axis) rotating at an arbitary ele- 
ctrical angular speed dO/dt. It is clear that sot 

is fixed in the stator and fixed in rotor 
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and honco rotates at an angular velocity Subscript 1 is 

used for stator quantities and subscript 2 for rotor quanti- 
ties. The time zero angular relationship botvj'een three set 
of axes can be- selected arbitarily; however, it is convenient 
to assume that at time zero ,3^ end ’d’ axes coincide. 

The transformation equations 'vvhich for the balanced systc 
can bo correlated to the angular relation of the axes shown in 


Fig. 1.1 < 

iro 
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as follows 
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In those equations the variable f can represent cither voltage, 
current or flux linkage. The equations are restricted in that 
the instantaneous angular displacement 9 of the arbitary 
roforcnco frame must bo continuous finite function. 

It has boon shovm in [4j that the equations of motor 
performance for squirrel cage induction motor for a stationary 
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refsranco frame are 
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(1.7) 


( 1 . 8 ) 


where m is the number of phases, P is the number of poles, 
is the leakage inductance of the rotor, L 22 ’the leakage 
inductance of the rotor referred to the stator, M is mutual 
inductance and ’p’ refers to ’d/dt'. In tiiOse 'T ’ is the 
instantaneous torque of the motor. 

In order to solve the above dq equations, the three phase 
input excitation currents are transformed into two phase dq 
currents. The transformation connecting 3 phase to 2 phase 
variables with assumed directioh of coordinates as shown in 
Fig. 1,1, for tho case of balanced stator current gives 


^dl ^al 

■^ql ^-^bl ^cl^ 


(1.9) 


( 1 . 10 ) 
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Axes of throe phase synxiotrical machino 
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CHAPTER 2 


STEADY STATE ANALYSIS OF THREE PM^\SE INDUCTION MOTOR 
FED BY AN IDEAL THREE PHASE CURRENT SOURCE 

IN’VERTER 


2.1 INTRODUCTION 

The three phase squirrel cage induction motor fed by an 
ideal square wave current has been analysed in this chapter. 

Fig. 2.1 sShows the waveforms of the line current in the motor 
produced by ideal current source inverter. The dq model of 
the induction motor, with axis attached to the stator (Sec.l.l)j 
is used for the analysis. In the analysis, the steady state 
solution has been obtained. This solution is for a constant 
invertor frequency and a constant rotor speed. The inverter 
frequency is adjustable by the inverter gating circuits. 

The three phase current is transformed into the d~q 
current using the equations (1.9) and (1.10). The transformed 
d-q currents, for all the six intervals of the invertor period 
are listed in Table 2.1. In Table 2.1, 

I = 2/3 I^^ (2.1) 

whore I^^ is the dc link current fed to the inverter. 

The waveform of the d-q currents is given in the Fig. 2.1. 
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The waveform of tho rotor current is obtained by the 
solution of equations (1.6). In tho case of current fed 
induction motor, the above equations are decoupled from eqn. 
(1.7) which correspond to tho phase voltage of tho motor. This 
is in contrast ivith tho voltage fed induction motor, in v/nich 
case all the four equations of the matrix equations (1.6) and 
(1.7) have to bo solved simultaneously. 

Since wo are solving for rotor currents at constant rotor 
sDOod, tho problem can be viewed as obtaining tho steady state 
solution of the differential equations with constant coeffi- 
cients. This has boon done using two well known techniques s 

(1) Storody state solution in time domain determined by 
obtaining homogeneous and particular solutions and using two 
boundary value conditions. 

(2) Steady state solution in frequency domain obtained by 
considering tho input as the summation of sinusoidal signals, 
through fourior series of stator current and calculation of tho 
transfer function for those inputs from tho equation (1.6). 


Equation (1.7) has been used in this chapter to obtain 
tho voltage waveform in tho d-q frame at tho stator terminals. 
Those equations have terms containing 'P^dl' *^^ql** 

'idi* and discontinuous (Fig. 2.1) the voltage wave- 

form is shown to have impulse. 
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Soction 2,2 deals vdth the calculation of the rotor 
currents, using the time domain analysis. In this the eases 
of the- stationary and the revolving rotors arc dealt with 
Separately. In ease of the stationary rotor, tho motor 
equations rela.ting tho dq stator and rotor currents, from 
eepjation (1.6) turn out to be first order and decoupled. The 
discontinuities in tho rotor currents duo to discontinuities 
in tho stator current can bo computed from equation (1.6). 


In ease of revolving rotor, equation (1.6) has two 
coupled equations for tho rotor dq currents, i^jj 2 and iqp* A 
single equation for iq2 obtained from equation 

(1.6) by elimination. However, tho resulting single variable 
equation is a second order equation which contains tho terms 


as ;5nd ’p^iq]_‘* Since tho stator. dq currents, i 

i^j_ arc discontinuous, these terms correspond to tho dcriva- 
tives of tho impulses. To avoid such terms, tho two now 
vaiiablos Xj and , have boon defined instead of the rotor 
current variables, 1^2 sind iq 2 » Those variables and X 2 
have boon shown to bo continuous and are proportional to tho 
flux linkage in tho machine. Those are sometimes referred to 
as pseudo rotor currents in tho literature [4]. 


dl 


and 


Tho solution of tho equation (1.6) in time domain 
requires two boundary conditions. Those have been obtained by 
noting tho fact that under steady state at each multiple of 
60"^ in ojt, to being the inverter frequency^ is similar. The 



st-’/tor mmf juraps in space by 60° due to the particular forin 
of i_,ii_ and i. (Fig. 2.1) vvhich have discontinuities . The 
stator ramf instead of rotating smoothly as for sine v/ave 
currents ju^aps by 60° in space for every T/6 period of the 
inverter p T being tho inverter time period. Thus, under 
steady state the rotor ramf must also move by 60° in space in 
each T/6 interval so that each interval is similar. The 
resulting boundary value relations based on above arguments 
are given in equations (2.5) and (2.5), 

Section 2 .2, the solution for the initial value of 
rotor currents has been obtained independently for each inter- 
val. From these it can bo seen the current solution obtained 
at the start of each interval is sane as that obtained for the 
end of the previous interval. Thus it is sufficient to solve 
for the initial value of rotor current in any one of the six 
intervals . 

Section 2.3 deals with the technique of frequency domain 
solution. For the rotor currents in this scheme- the expre- 
ssions for tho rotor harmonic currents have been obtained in 
terms of the harmonic components of the stator current. The 
cases of stationary and rotating rotor have been dealt with 
in ono general expression. 

In Section 2.4 , the voltage expressions in dq frame at the 
stator terminals have been obtained in time domain. The results 
of Section 2.2 and matrix equation (1.7) have been used for the 
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c ,-Lcuxr': 'Cion 


Foi 


verification of tho solutions through 


to uh tlio 52clini:u:;s for rotor currents, tivo computer proprr.a 
h u''G boon developed. Tne first propreei computes the tifac 
do.Te.in solucion of Lhs rotor current and then obtains the 
hercuionic coiaponcnts of this solution. The second pjrojr.ims 
directly cououtes the he.r ionic coiaponents of the rotor 
cui'j'ent by r.l'.o frequency doaiein tcchnicaie (Section 2.3). In 
Section _2.5 the ceso of r> p'lrticulnr induction notor hos been 
studied. Usinu; tl'io p ’r.'i’ va ters of this induction iTiOtou', it has 
b-'en si'own tbrt results obtained through both programs are in 


:lo; 


iroe.oorrc, 


This verifies these x'esults. 


2.2 Ai'.IALYSIS Or INQUCTIOM MOTOR FOR ROTOR CURi-^EMTS THROUG'I 
Ti; i-i DOM/JN 


In this section, tho steady state solutions through tho 
tirao do;i''in technique, Section 2.1, for the rotor currents 
have -aoen obtained. For this the motor perf oraance aquation 


(1.6) has boon solved. The c^ses of stationary and the 
rotating ro to.r hove boon c’ealt with separately. The equation 
(1.6) lias boon solved for the first throe intervals (I, II, III) 
o < (.ot _< Tc . The stator current waveform (Fig. 2.1) is 
symiiiatric about cot ■= n. Therefore, solutions for intervals, 
IV, V and VI, it < oot < 2-ic,aro negative of the solutions of tho 
intervals I, II and III respectively. 
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For tho analysis, each interval has been associated v/ith 
the discontinuity at the start of the interval, '.7hile solving 
equation (1.6), each interval has boon divided into two sub- 
intervals. One corresponds to the point of discontinuity at th 
start of tho interval and the other to the reinaining period of 
that interval. 

The following syiabols have been used for the analysis. 


Ido = Id2(h = 0") ; = 0-) 


Id2 " ■^d2' ‘"i 


= 0-^) ,• IJ,2 = lq2(h = O'") 


= TJ ; 1^2 = = Tc) 


■d2 " ■^d2'''l - ‘o 


( 2 . 2 ) 


where T- is the duration of each interval, i.e. (T = n:/3(o) . 

The tirae is the time within the interval whore i denotes 

the interval number i.e. i = 1,2,3,, Thus tj^,t 2 and t^ are 
defined for intervals I, II and III respectively and in these 
we have ( tj^ = t) , (tq = (t“n;/3(o)) and (t^ = ( t“-2'Ji/3L0 ) ) . As 
already iiientioned that each interval is associated with the 
discontinuity of stator current at the start of that interval. 
Thus ( t^ =0*) refers to time just before the discontinuity 
has been considered in interval ’i’ and ( t^ - O"^) corresponds 
to the tii.ae just after the discontinuity has been accounted for. 


2.2.1 Case of Stationary Rotor 

Frora equation (1.6), the equations relating tho dq stator 
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snc; rotor currents for the stationary rotor 

i'AFi^l + (^2'^^22 ^^^c12 ~ ^ 

iApiql + (Y'2+L22p)iq2 “ ^ 

In this sub-sGction we are first obtaining the solution of iq 2 
and i^i^ from above aquations with the arbitary initial values 
of and Then, the fact that steady state dq currents 

at tlio G'ld of each invertor period make the airgap ramf v;ave 
turn %/3 electrical radians, is used to compute and 

The values of dq currents at the end of the inverter 
period are related to the initial values by the following 
equations [4], [ll] 

Id2 = (Ido - V3 Iqo)/2 (2.5) 

Iq2 = (1,0 '/■3do)/2 (2.6) 

where 1^*2 snd defined in equation (2,2). 

In Section 2. 2. 1.1, the rotor currents expressions for 
interval I have been computed. Sections 2, 2. 1.2 and 2,2, 1,3 
give the rotor current expressions for intervals II and III 
respectively, v/hich can be obtained by analysis as done in 
2, 2. 1.1, Sec, 2. 2. 1,4 shows that it is sufficient to calculate 
Ij^Q and IqQ for one interval only. 


= 0) are, 

(2.3) 

(2.4) 
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2. 2, 1,1 Interval I {o < cot < %/3) 

As tj^ — t in this interval), the time is denoted as t 
i n th is interval. 

CaJLcula;tipn^ of currents _at t_^ =^0 

From Table 2.1j at t = 0> 

Ai^^l(t=o) = i^^(t=0"^)“i^j_(t=0“) =0 (2,7) 

Ai^jj_(t=o) = = “fSI (2.8) 

JL, 

Integrating equation (2,3) from t = 0 to t = O' , 

O"'’ 0'*' 0"^ 

M / (pi^j,)dt + ^22 ° 

t=0" t=0“ t^O”* 

Since i ^2 finite, this gives 

So frora (2.7) 

Ibiicii (t=0) ^22^1^2^'^=°) = ° 

So frora equations (2.2), (2,7) and (2.9) 


^d2 ^do 


(2.10) 



2.9 


Similarly fro;*i (2.4) using {2, A) ond (2.3) 


i^^(t=0) = 




-22 


•ql 


S O y 



+ f3 Kj_ 


X'lh'd'CQ 


‘"1 


ll.-, 

L^9 


-alculntion of currents for 


F3:or Table 2.1, during 


t >_ p 

this period 


idl = --31/2 

iql = ■■'.'■31/2 

Solving (2.3) and (2.d) using (2.14) and (2.15) 

“ ^d2 e-P(“a‘t) 
i^^o(t) = exp (“at) 


V'/here 


a 


rp/L 


22 


(2.11) 


( 2 . 12 ) 


(2.13) 


(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 


and and 


are as defined in equations (2,2). 



‘5.°". .°A ido. 


SuUstitu cincj (t = Tt/Sto) in equations (2.16) and (2.17) 
ad usincj (2.10) and (2.12) we have 


T ? ? 

■^d2 


Ido exp(-aT^) 


(2.19) 


I"o '-= (I ^+’>r 3 K,) Gxp(™aT^) 
q.i ' qo ‘ 1 ■ ^ c' 


( 2 . 20 ) 


Solving for and frora equations (2.5);, (2.6), (2.19) 


nod (2.20) 


qo 


'CIO 


V ? 1- ’ 

3 ‘^1 ^"2 

d' ' 


( 2 . 21 ) 


qo 


-^3 kA 

[exp(-”aT^) 


1 / 2 ) 


( 2 . 22 ) 


v'/'iiero 


4 


exp(-aT„) 


(2.23) 


exp ( -2 aT ) -exp ( --aT_ ) +1 


(2.24) 


2. 2. 1.2 Interval II (i.e. | < ut < or o < t 2 < §) 


The rotor current results have been obtained for this 


interval preceding similar to Sec, 2. 2. 1.1, 



2„11 


At = 0 


■ * = T 

'd2 "^GO 


P 'c ’ 
2 1 


“if 

I'., = I ^ -Sv k* 

q2. qo 2 1 


For ±2 > 0 


(t^) =1^2 e-pC^a'to) 


d2' "2- 


(2.25) 


iq2('*^2) = Iq2 ©xpC^-atq) 


( 2 . 26 ) 


ana 


ir ? 


K 


^do 


= I - [Gxp(~aT^)l 


(2.27) 


I 


qo 


Y3 I<] exp(-aT )^ 

[1 iA] 

^4 


(2.28) 


2 % 


2.1.3 Interval 111 (i.e. < t < tt or o <_ t^ < 


Similarly, for this interval. 


at to = 0 


^d2 ^do “ 2 


I > = I ™ 

q2 qo 2 1 


(2.29) 
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for to > 0 

O 


= ^d2 


(2.30) 




Iq2 exp(“C»t3) 


(2.31) 


anc 


’do 


qo 


R ' , ’ 
3 1-2 


C, 


:• [exp(“aT^)"l)] 


■ r . 


. K’ 

I [exp(-aT^)+l] 


(2.32) 


2. 2, 1.4 Discussion 

The rotor current expressions for intervals IVjV and VI 

can be observed directly fror.i those in S. e^-M-en - c 1 ,11 and III 

respectively because of syrivaetry property, Ue also observe 

t'nat for this case of stationary rotor, the closed form 

expressions for 1^.., and I ^ are obtained. In the above 

c V Q q o 

sections, the analytical expressions of I^^ and 1^^ have been 
computed independently for each interval. The expressions 
of the final values, at the end of each interval, can be 
obtained from equations (2.5) and (2.6) using the expressions 


of I 


do 


and I 


qo 


obtained for that interval. These values are 


suipjriarised in Table 2.2 


It is seen frora Tables 2.2 that the rotor cur.rents at 
the start of each interval i.e. at (t^ =0) is same as that 
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o':/‘c?i;'ioc'; at the earl of previous interval i,e. at ~ 

Thus "/c need to solve for and I ^ for one interval only. 

GO Q ^ 

Case of Rotating Rotor 

Fi'Oii equation (1,6), the equations relating the dq stator 
cO. n rotor cui:-"onts r/ith '.o, as the angular velocity of the 
■ /:or in oioctricai radians per second are 


.’lU 


r "-dl 




2 "22^'" ^ ^c2 ^22 r cl2 


(2.33) 


I'jI 


+ to., 1 J_ + ^^2'^'^22P^^d2 “ 


0 


(2.34) 


These Q-^v hions ore coupled in contrast with tho stationary 
ro !:or case, equr'i tions (2,3) and (2.4) ’.-/here these are decoupled, 
I’urth^'r pi(j] 3'id piq] aj-'e impulse functions os and i^j^ are 

dlscon ti'iuous , The solucion due to iiOpulse functions has been 
d -n -0 fallows 


To obtain the changes in rotor currents at any tiiao instant 

’t-./ due to change in stator current va'e integrate equation 

(2.33) from t = t. to t = t” . This gives 

o o 


o 


4- 

t; 

•4 




''V ^■dl'^’^ + 


t=t 


o 


t 

o 

r> 

J 

t = t‘^ 


t 


i^gdt + L22 


/ p i.oCit 


t=:t 


■q2 




L/^rj J Cl^ ^ 


(2.37) 


o 



2 . l-i- 


iirc-nts ero rinite, 


It = 0 wh: 


^ 2 ’ ^'-1 


0.-, or ci. 


(2.30) 


, from (2«37), (2.30) 


,[i ,(t-f^) 

C|1 0 




C;2 O 


=f:)] = 0 


(2.39) 




f ' -}- 'N 

■ql^ 


(2.40) 


whore (t=t^) corresponds to any arbitary time instant and ’ 
corresponds to a jump in tbo current variable. Similarly, 
from equation (2.34) 




4;- 4i(t=tp 


(2.41) 


Thus the effect of tho discontinuity in i^^^ and i^j^ at 
tj^ = 0 is a corresponding discontinuity in i ^2 ^q2* 

by equations (2,41) and (2.40) respectively. This is equivalent 
to stating tliat tho response due to 'pi^ 2 _’ find ’piq^' impulses 
in equations (2,33) and (2.34) can be looked upon as sudden 


changes in the initial conditions of i^^ and iq2' 


Therefore, 


■t is possible to solve (2,33) and (2.34) as a set of siraul- 


For the interval t^^ > o, because i^j^ and 


tanoous equations. 
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"ql 
a s 


are constants (Fig. 2.1) and pi(j 2_9 P^ql zeros. Thus 

olution similar to Sec. 2.2 of stationary rotor can be 


9 


obtained . 


Hov/ever, equations (2,33) and (2.34) can be solved more 
conveniently if a change of variable is done from i^,^ and i^^ 
to and respectively. These variables, and are 
defined as, 


''^1 “ ^q2 ^ql 


(2.42) 


-2 = ^d2 "■ 


,iL i 
1 / 2 2 ^ ^ 


(2.43) 


It is shown in the following that these now 
variables and X 2 are continuous functions of time. 

The change in X^ at a arbitory point (t = t^) can be 
obtained from equation (2.42) as 


>M(t=t^) =-iq2('t='to)+ T™’^ ir,i(‘t=t^) 


■‘22 


ql' 


(2.44) 


Using (2.40) in (2.44) wo get 


Zi = 0 (2.45) 

Similarly from (2.43) and (2.41) 

A-X2(t=t^) = 0 (2.46) 

These new variables X^ and X 2 referred to as pseudo- rotor 

currents [4] and arc proportional to the rotor flux in the 
machine . 
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In equations (2.33) and (2.34), eliminating feear ’P^d2* 
and ’piq 2 * ’^y using (2.42) and (2.43), we obtain 


P-^l = ^r^^2 " ®iq2 


pX2 = “ 31^2 


wneue 


- 2'''^22 


(2.47) 

(2.48) 

(2.49) 


Substituting for i^q and i^^ in equation (2.48) and (2.49) from 
(2.4^) and (2.^jfe) gives 


“■I ~ r"^2 


^""1 Lr 


(2.50) 


pX^ =-^to X, 


®"^2 L22 ^^dl 


(2.51) 


50 7 

Dif forontiating (2.4^) and substituting for pXo from (2.44) and 


similarly differentiating (2.51) and substituting for pXj^ from 
equation (2. 48), gives 


2, . 2^.„ 


where 


= qaiql+q^r^dl+’^lPqi 

(2.52) 

= Kj^alai-Ki“rqi+4Pqi 

(2.53) 

— A^ia /^2 g 

(2,53a) 

and (2.53) are the final equations 

rotor current expressions. 

which are 
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2. 2, 2,1 Interval - I ( o < cot < %/3) 

.-.t cot = 0, froro continuous property of X-, 5 X 2 

= 0 ? ^ X 2 = 0 (2,54) 

for 0 < cot < 7 i/ 3 , fro;i Table 2.1 


i 

11 

1 — 1 

x> 

•H 

-31/2 


^ql - ' 

-f3l/2 

(2.55) 

Frora ( 2 , 52)5 for 

this interval 



Q r) K 1 B'\f 3 1 K T CO 3 JL 

(p + 2 ap+a +co^) X^ = — . .. 


The solution of above is 


Xj(t) = exp(»at) ^ol (2.56) 


where 






'2 ‘7 '2 
(a +cop 


(2.57) 


a no. 


C’ , D' are constants, 
ol' ol 


'■10 

] 

= ^^ 1 ! 
1 

1 

t .=0 

J 

''20 

— 74 

- --'2 

1 

t . =0 

1 


• 

''10 

cl X - 
== dt ' 

L 

0 

9 

0 

Mo 

6 X 2 1 

“ eft '1 t^=c 

(2,58) 

X” 

"'1 

i 

- "Ml 

JL 1 


0 

9 

Y *'* 
^'2 

— Xrs 

t,= M 
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where i corresponds to interval. Hero i = 1 and tj^ = t. 
Using above in (2,56) gives 


°ol “ 'So * *^4 

(2.59) 


(2.60) 


Solving (2.50) at t = o, using (2,55) and (2.58) 

^10 = ^‘'r^20 " "^'"lO (2.61) 

Similarly solving (2,53) for X 2 gives 

X 2 (t) = oxp(«“ 0 t) [C ^2 “r'*'^^o2 (2.62) 


where 


= 


V3 

5 


K3_I(w^- 
( a2+u)p 


a’fS) 


^o2 “ ^20 " ^6 
^o2 = !^^20‘‘'^°o2^/“r 


X, 


20 


■to X , 


r 'UO 


“ ®^20 "■ "S ^1^ 


Computing X^^ at t = from (2.62), using (2.58) 


V 


t=T 


= '^il ^3 ^ °ol ^4 


K 


■4 


vdiore 


Yo = exp(”aT ) sin to T 
O C X c 

= oxp(~3T^) cos to^T^. 


(2.63) 

(2.64) 

(2.65) 

( 2 . 66 ) 

(2.67) 

(2.68) 


(2.69) 
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■Substituting for in (2,67) from equations (2.59) 

and (2.6C). Using (2.51), we have 




where 


-1 = VlO ^3'^20 ^3 


Z’ = C -/3 - 


= [- K^I + K^aj/o)^ 


(2.70) 

(2.71) 

(2.72) 


Similarly computing at t = and solving as above we get 

(2.73) 


y !? _ _Y Y -r V Y 4- 7 ’ 

■ ^2 “ ^3"^ 10 ^4" 20 4 


vdiero 


where 


Z 


4 


r ^ 
^6 


^ 6^3 ” V 4 ^6 


[» I Kj_I - aK^]/a), 


v6J/-r 

.■/riting equations (2.70) and (2.73) in matrix form 


(2.74) 


(2.75) 


1 1 

1 

— 

— 


— ~1 


— 

! ? ’ 

^2 ■ [ 

j 

1 

i 

“Y3 




71 

^4 

1 YU ! 

! 1 1 


I'S 

^4 


. ^10 


71 

^3 

j i 


1 



J 




(2.76) 


From the boundary condition property, the relation between 
the initial and final values of current of an' invertor 


interval, 


X-> 


1/2 

f3/2 


~f3/2 

1/2 




X, 


20 


V 

A 


10 


(2.77) 
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From equations (2,76) and (2.77), solving for 


^20 


- Y4) 

-(^3 . y3)- 

-1 

N 

1 

^10 


("'' I - 

- Y4) 


7 ’ 

^3 


X 2 Q and known from (2.78) after calculation of constaii 

in the order K^, Y 3 , Y^, C^, and Z^. 

The pseudo currents X^^ (t) and X 2 (t) are now known from 
equations (2.56) and (2,62), The actual rotor currents i ^2 
iq 2 can bo obtained from equations (2.44) and (2.45)* 

2. 2. 2. 2 Interval II (i.e. | < wt < or o < t 2 < 3 ) 

Solving as interval I we obtain 


Xj^(t 2 ) = exp(“at 2 ) [cos 

sin “f" ^q 3 

(2.79) 

where 



. ^ T3 K, la 

7 

(a +CD^) 


(2.80) 

^o3 = ^10 ^ 


(2.81) 

'^o3 ^ "^10 ®°o3 


(2.82) 

and 




^ 2 ^^ 2 ^ = exp(“at 2 ) [Co_^ sin ^^p‘C 2 '*'Do 4 


(2.83) 


where 

K 


10 


f3I Kj_to^ 
(a + 0 )^) 


(2.84) 
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%4 ^20 " ^10 


*^04 ~ sDp,^]/co 


"20 


04- 


(2.35) 

( 2 . 86 ) 


10 > -^10’ ■^"10 ^20 ^°^^®spond xo 


It should be noted that X^q, X, X,^ and 
definition of equation (2.58) with i = 2. 

For initial condition of this interval we have the 
relation 


^"20 

! 


(-1 - Vp -(\’| - Yj) 

-1 

z, 1 

8 I 

% 

o 

t — 1 
>< 

' * j 

_ (''2 - Y3) (| - Y4) 


I z- 
L 7 J 


(2.87) 


where 



= 

Ci2 ^3 " K10Y4 

^10 


(2.88) 

Zij 


Cll ^^^3 + ‘<7 ^4 

" K 7 


(2.89) 

^12 


“%0 



(2.90) 

^11 


[~f 3 Kj_I -y K^a] 

/“r 


(2.91) 

X 

0 

and 

X^Q are solved, 

the pseudo 

currents X 2 _(t) 

and 


X^it) are defining over the entire interval %/3 to 2%/3, The 
actual rotor currents can be found out using equations (2.44) 
and (2.45) , 
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'TT 

2, 2. 2. 3 Interval III (i.e. 3 ” 5. ^ ^ 0 ;< < ■j) 

Solving as for interval I we obtain for this interval. 
X^(t) = exp( -at^) [cos ' sin t^ ■+ cos + ^^^3 (2.92 


X2(t) 

= exp(~at 3 ) [Co^ sin w^t^ + cos w^t^] + 

(2.9f 

where 




= t - 1%fZ 

(2.94) 

^13 

fSKj^I (f3w^-a) 

( a +0)^) 

(2.95) 

''le 

V3K,I (to +f3a) 

_ fa 

(2.96) 

Do5 

“ ^10 “ ^13 

(2.97) 

°05 

= [X^Q + aD^^]/(A)^ 

(2.98) 

Do6 

= X 20 ” 

(2.99) 

*^06 

= [X 20 + 

(2.ior ' 


For the initial condition, the following relation is obtained 


^20 


(j -* Y4) 

- ¥3)' 

-1 

" 2 : 12 " 

^10 • 

j-s. 

i 

- V3) 

(|- - Y4) 


2^ii 


( 2 . 101 ) 



2.23 


where 


Zi2 


= 48 ^3 “ 46 4 1^16 


^li ~ 4? ^3 “ 43 4 43 


where 


, ^3Kll 


r 


OT/ T 

C«8 = [ ^24 - Ki^aJ/o)^ 


( 2 . 102 ) 

(2.103) 

(2.104) 

(2.105) 


Thus, the actual rotor currents have been found out for this 
interval. In rotating rotor case it can be seen that analytical 
expressions are complicated and closed form expressions are not 
obtained for dq rotor currents. A computer program has been 
developed to calculate initial and final values of currents 
for these intervals in rotating rotor case. The boundary condi- 
tions have been shown- to be matching for the case of induction 
motor in Section 2,5. 


2.3 ANALYSIS OF INDUCTION MOTOR FOR ROTOR CURRENTS THROUGH 
FREQUENCY DOMAIN 

The expressions for different harmonics of the rotor 
currents have been computed in this section. The cases of 
stationary and rotating rotors have been dealt with in one 
general expression. 
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The analysis proceeds in two s teps.F irs tly , is the 
calculation of the harmonic conponents of the stator current. 
Then equations (2,52) and (2,53) are used to compute the 
harmonics of pseudo rotor currents and With this 

result, using equations (2,42) and (2.43), the harmonics of 
rotor currents can be obtained. 

For calculation of the harmonic components of stator 
current, i^j^ and i^j^ (Fig. 2,2) are resolved into Fourier 
series as 

CO 

i ,, = S a cos nut (2,106) 

n=l "^d 


i , = S b sin nut 
n=l 


where 


for all even values of n 


a =0 


bn =0 

Cl 


and for all odd values of n 


61 . nil 

1 = sin 

n^ n-rt 3 


K = . 4m ™ 

n nit 3 

q 


(2.107) 


(2.108) 

(2.109) 


( 2 . 110 ) 


( 2 . 111 ) 
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It can be seen froui (2.110) and (2.111) that 
triplensi.e. n equal to a multiple of 3, a 


Thus 


for all 

and b = 0, 
n 


f or n 


a 


n 


d 


b 


n 


q 


1,5,7,11,13,17, ... 

(2.112) 

- -H- 

n7i 3 

(2.113) 

„ K3 j 
rm 

(2.114) 


and all other harmonics are absent. 


In the remainder of this section only o dd and non triolent 
h armonic s are considered from (2,106) and (2.107) 


. 61 . nit ^ 

^d " ^ “ ni" 

1 n 


i = S “ sin nwt 

Qt riTi 

u n 


V/riting for mth harmonic of i. and i 

d^ Ql 


qi = ^md + 0!^j) 

m 


cos(™t + 


So, from (2.115) and (2.116) 


(2.115) 


(2.116) 


(2.117) 

(2.118) 
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Imdl 

= l^mqt = llml 

(2.119) 

^'md 

= 0 

(2.120) 

a 

raq 

= ~ %/2 for m = 1,7,13,19, ... 

(2.121) 

a 

mq 

= + %/2 for m = 5,11,17,23,... 

(2.122) 


VVe define a variable *p' as 


1 if m = 1,7,13,19, ... 

P = (2.123) 

0 if m 5,11,17,23,... 

then from (2.117) and (2.118) 


^dlm = 

^qlm ^/2) 


(2.124) 

(2.125) 


wh e re 

W- = ®nd ( 2 - 126 ) 

Solving (2.52) for rath harmonic of X-j^ denoted by equation 

becomes 


( p^+2ap+a'"+oj^) 


K, ai T -I- K,(jo i ,1 + K, pi , 

1 qlm 1 r dim 1^ qlm 


This equation can bo written in phasor forra as 



where 


At 

r, 2 2^ 2^ 2.2^,^ ^2^l/2 

= [("m CO +a "Ho^) +(2atom) J 

(2.128) 


= tan~^ [2amco/(a^-Ko^“ra^co^) ] 

(2.129) 

^qlm 

= ^ /ic:l)P I 

(2.130) 

^dlm 

= ii! /o° 

n ~ ■ 

(2.131) 


and corresponds to the phasor of from (2.127), 

(2.130) and (2.131) 

K,a K,oj 

Iqln> '-X" ^dlm + 

mwK, ^ 

+ (-^^ ~ a^) (2.132) 

-- > 

This can be visualised as if X,^, is obtained from summation 
of three phasors which can be obtained from the aod 

by passing these through apprqoTi ate amplifier, as shown in 
Fig, 2.2. 

To obtain ^^^^^^(t), taking the real part of equation 
(2.132), we obtain 
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K al 

Xij^(t) = — • cos(mwt--aT,+(-l)P -11/2)+ — cos(iniji)t~aT.) + 


Ap 


A. 


mu)K, 


"A^ + I + (-1) I - “t’ 


This gives 


X 


K3_aK)^ 


Im 


K,! p 

siri(ma)t-ap)+ — (a)^-(-.l)^ mw)cos(niwt-ctj.) 

(2.133) 


Similarly we can obtain mth harmonic of X 2 j defined as X 2 j^ 
from equation (2.53) as 


X 


K^a(-l)P 


'2m 


A-t 


m 


cos(m(iot-a^)- 


Ibin 

At 




sin(mmt-a^) 


(2,134) 


It is noted from equations (2,123) and (2,134) that magnitudes 
of X^j^ and X 2 jj^ are equal and they have a phase difference of 

Tt/2, 


Thus the harmonic components of the pseudo rotor currents 
and X 2 are computed. Using equations (2,44) and (2.45) the 
harmonic components of the rotor current can be obtained. 

Using the above procedure, the harmonics of the pseudo- 
rotor currents have been conputed for a case of particular 
motor, in Sec. 2.5, 
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2,4 VOLTAGE V/AVEFOPJ';! AT THE STATOR THR'IIUALS 

In this section, tho voltage waveform, in the dq frame, 
ot the stator terminals has been computed from the equation 
(1.7) , 


Vdi (r^+Ld3_p) ^dl ^'p^d2 


\il = ^^l-"4lP^ Si ^ 


(2.135) 

(2.136) 


The expressions for i^^ and 1^2 obtained in Sec. 2.2 have been 
used in the above equations to compute the voltage y/aveform. 

The casGS of stationary and rotating rotor have been dealt v/ith 
separately . 

2.4.1 For the Stationary Rotor 

For the stationary rotor, substituting for pi^o from equa- 
tion (2,3) into equation (2.115) we have 


\r 


dl ~ P^dl 1^2 ^d2^^ 


‘■2 . 




Using (a = from equation (2,18), in above, 

\ll ^l^dl“^'^® ^dl ^^11 “ p^dl 

Similarly, from (2.4) and (2.136) 


Si " n^ql‘'^'^®^q2^*^^ll” 122^ 


(2.137) 


(2.138) 


It can be noted, thiat for ideal square wave current, pi^ij and 
pi^^-j are ' j, except at the end points corresponding to 
cot = 0, it/S, 2 tx/3 


y • o i 
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At t 


!" OSG point, has impulse. From equation (2.117) 


if. is observed tha !; at these points of cascontinui by of i^i-,,'/ 


1^ ''dl 


S (t^=o) 


''■'here 


nagni 

tudo 


'^1 


(2.139) 

Ai = 


(2.140) 


% 


tl'.ere 

is an iiioulse in i 

:h e voltage via v e f o rm 


because of the leakage inductance. 

It is also observed from equation (2.137) that there a 
j jmp in also apart fro.n i'.rpulso because of discontinuity of 

idi ?-od idp* hy 'a ' we refer to the change at the point of 

consiciGrationi, from (2.137) vie have 


AVdi(ti=o) - I'lA id]_( tj, =o)~Ha/A id 2 ('ti = o) (2.141) 

Su’ts titu ting forZ\idp ~ froia (2.9) into (2.141) we have 
AVdi(td=o) = idi(%=o) (2.142) 


whejee 


A 


.^2 




(2.143) 


Because (aar 2 ) and honco this jump is because of rotor and 
stator resistances v^/hich should bo so as there is a suddon 
change of cu;:rent through those resistances. 
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Similarly s, for wg have from (2.138) anc’ (2.11) 


>' ( t . =0 ) 


v^i = *1 Ayf y=°> 


\ V . ( t . =o ) = Ao A i , ( t . =0 ) 
■* qi 1 ^ ^ — ql^ 1 ' 


(2.1/14) 
( 2 . 145 ) 


During th? time ( > o) Vs/e note that i^^^^ and are 

constants and hence equations (2,137) and (2.133) reduce to 

^dl ^'l^dl ” ^ ° (2.146) 

V “ ^4ql “ h ^ ° (2.147) 

Thus, hno^'dnq stator and rotor currents expressions, the 
voltages and v^j^ caii be computed the impulse and jump 

magnitudes can be computed from the information of discontinui- 
ties of i^j^ and ijjj_. Table 2,3 gives the reference equation 
nuiabers which are are to be usod to obtain voltage expression 
by use of equations (2.142), (2.144), (2.146) and (2.147). 

The table is only for interval I, II and III because expressions 
for intervals IV 5 V and 71 can be obtained from these because of 
symmetry about (ojt = 71 :). 

2.4.2 For the Rotating Rotor' 

From equation (2.33), we have 


pi 


q2 


“r^d2 


^^-q2 1^2 ^r^dl " P^ql 


^22 
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Substituting for 'pi^,^’ in (2,136) from this equation. 




(2.148) 


where 'Aj_’ is given by the equation (2.140). 
Similarlys, from (2.34) and (2.135) we have 


"^dl ” ^l^dl"^''^®^d2‘“^'*'^^r^L”2 ^ql'^^q2^'‘'^l^^dl 


(2.149) 


From (2.148) it can be said that at the end points of the 
intervals (wt = o, 'ir/3, 2 tc/ 3). has an impulse and a jump 

whose expressions are given as 


6 (t^=o) 




(2.150) 


^ (2.151) 


Substituting for i ^9 in terms of iq 3 _ '^^d 2 terms 

of ds i(^ 2 _ from equations (2,9) and ( 2 . 11 ) respectively, in 
(2.151) we have 


= *2 A lqi(ti=o) 


(2.152) 


where A 9 is ’given by equation (2,143), 


Similarly from (2,149) using (2,9) and ( 2 . 11 ) 


AV,,(t. = o) = A, ^ii,,,(t. = o) 


(2.153) 


(2.154) 
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Thus, the behaviour of and at the points of 

discontinuity of i^^^ or i^^ can be computed from the equations 
(2.150), (2.152), (2.153) and (2.154). 


For the remaining period of time, we note that and i^^ 

are constants, Fig. 2,2. So, for t_^ > o, equation (2.148), 
(2.149) reduce to 


V 


ql = ^l^ql ■ ^dl + 


'^dl ^l^dl " ‘''®^d2 


From equations (2,42) and (2.'^!3) 


(2.155) 

(2.156) 


^q2 ~ ^1 


1 - 9 ^ ql 


(2.157) 


■d2 " ^ ^dl 


(2.158) 


Substituting for i ^2 equations (2.157) and (2.158) 

in equations (2.155) and (2.156) we have 


'^dl = ^ bl ' ^ - ‘■'“r 


(2.159) 




(2.160) 


Since X-j^ and Xp expressions have been computed in Sec, 2.2.2, 
the expressions for and ^^n be computed. 



Table 2.4 lists the reference to equations numbers of 


and during intervals 1,11 and III v/hich are to be 
used in equations (2.159) and (2,160) to compute and ^qi* 

The table also lists the reference to changes in currents at 
t^ = 0, which can be used in equations (2.150), (2.152), 
(2.153) and (2.154) to compute the magnitudes of impulses and 
jump at (t:i = o) for the waveform of and 


ql* 


2.5 CALCULATIONS FOR A IQiO'm I-iACHINE PARAIITaRS 

For the calculation, a slip ring induction motor 
considered in [l2] has boon taken. The parameters of the 
machine are given in Appendix A, 

For this machine tv^o computer programs are developed in 
DEC Systom 10 computer in FORTRAN. One uses the results of 
Section 2.2 to compute the sampled rotor current waveform 
in time dom.ain. The other program computes the various 
harmonic components of the pseudo rotor current using 
Section 2.3 approach. The listings of those programs 
are given in Appendix B. These programs have inverter 
input dc current, inverter frequency and rotor speed as 
parameters . 

Table 2.5 gives the values of currents at the start 
and end of intervals I, II and III using results of Sec, 2, 2, 2 
for rotating rotor case. It is seen from this table that the 


rotor currents at the start of eaett jinteryal i.e. at (t.=0) 
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is same as that obtained at the end of the previous interval 
at - '^c^ * Thus and are to be computed only 

for one interval. 

Table 2.6 gives the harmonic components of the time 
solution of pseudo rotor current obtained through first 
program. The harmonic components which are obtained by 
second program using Sec, 2,3 have also been listed. Both 
those components are seen to similar and this varifies that 
both the programs give same solution. 
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Table 2.1 

Throe phase and two phase stator currents 


Interval/ 

conducting 


o < wt < n/3 


3 phase currents 


D~Q curr-G^vtsr 


II 


a V . 271 

e Wt \ 


III 

< 0}t < 71 


■al 

'bl 

cl 


■al 


■cl 


'al 

•bl 

■cl 


‘^dc 

0 

^dc 


'dc 

"dc 


= I 
- •»! 


dc 

dc 


0 


'dl 


--3I/2 

iql = “V3I/2 
(S, i|^j^(wt=0) *= 0 

^qi(tot=o) 

^dl “ ^ 
iqi = -r3i ^ 

•^^dl(i*>t=jr/3)‘"^^/^ 


icJi = 31/2 
iql = ~r3l/2 
Aidl(a}t=2it/3)'*^^'^^ 
^ ^ql(wt=27f/3)*^^^/^ 


IV 


Tu < a)t < 


471 


V 

47C X ^ Sti 
3“ < wt < ^ 


VI 


571 


< wt < 271 


■al 


1, 


bl 

■cl 


i_-> = 


■al 

■bl 

'cl 


'al 

'bl 

cl 


^dc 

id = 31/2 

0 

iql = T31/2 

‘^dc 

A^dl(tA)t=n) “ ^ 
^ ^ql(a}t=it) ”* ^ ■ 

0 

idi = ° 

^dc 

iql =731 

“dc 

^^dl(a)t=4ii/3)~ 

^ql(wt=47t/3)* 

■Idc 

idi = -31/2 

^dc 

y = r3l/2 

0 

A^dl (wt=57t/3)“' 
A ^ql(oi)t=57c/3)~' 


“31/2 



Initial and final valuos of D-Q rotor currents for stationary motor 
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Reference for the computation of stator voltage expressions 
in the case of stationary rotor 
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Reference for the computation of voltage expressions of 
rotating rotor case 


i 

5 

j 

ff 

f 

r-H j 

j 

T5 cro 

iH rH iH 1 

•H 'H j 

♦ • • f 

. » 1 

CN 04 04 1 

-H i 

O 0? o f 


**-4 r~l rH 1 

1 T3 ' 

x> f 

1 o C -p t 

<13 ft3 <T3 1 

1 ‘H ftJ (U j 
t 

1 

1 ! 

H H H I 

? 

J 

1 

f 

^ ^ ^ 1 


CM CO CO ! 


vO OD O 1 


t * « I 

CM CM CM 


"S— > S 

u 

* 

o 

• • • 1 

«+-» ; 

c c c 

a 1 

or cr cr ; 

1 

! ! 

1 

' 1 

O CD o 

i 

1 

i ^ ^ ^ 1 

1 

nO O CM 1 

5 

tT) Ov 1 

1 

» • • 1 

1 X 

CM CM CM 1 

1 

V.— ^ Nw^ j} 

1 U 

• • 1 

1 

c c d 

1 ^ 

cr cr cr 

1 

o o o 


1 

1 

1 

rH r-4 

... 

TJ tr 

2 

2 

2 

•H 



O 0> <D ! 

U *0 

rH rH fH I 

O C 

Xi io Xi 


CD <0 ro ! 


H H H ! 

I 

1— 1 

! 

ftS 


> 


u 

M 

CD 

M M M 


tHl IH 

c 


M 




TABLE ?,5 


I.H.V,Bff; 


=*'+EOK XI 
CHTTIAT. j 

VALUE ^ 


FINAL 

VALUE 


FOP 10C = 6.AMP,|ljyv 


2.96724 

■0.13596 

■3,10320 


■-kSim : 


FOR lt)C=6.fHP,nNV. FRE0,?28.HZ i 


“T:5n^r- 

III -'0,9299« 

for I0C=B.A«e ;IM. 

II «;5Jo^o 

m -0:88^3 

Hi 4 :imi 


o“44|9{"‘ ■ 

‘•''gOjJlOjHI R 

«0, 87940 
-■J, 39201 ■ 

-} .36982 

-1.59854 


**fnH K2** 

initial 

VALUE 

RPM=500. 

-1,87004 

-3.50471 

-1,63467 

RPM5300, 

•0,2807B 

-1,32993 

-1.04915 

RPMtsSOO. ' 
■0.21476 - 

•1,31291 
'1,09815 

RPM=300. 

■0.66n3 

'1,71524 

'1,05352 


final 

VALUE 


-3,50436 
— 1 ,64064 
1,87002 


•1,32837 
■1 ,05037 
0,29079 


■1.31125 

■1,09935 

0.21476 


■1.71421 

•1,05581 

0.66172 
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91.1185 
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90,9977 
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CHAPTER 3 

STEADY STATE ANALYSIS OF THREE PH,\SE INDUCTION MOTOR 
FED BY A CURRENT SOURCE INVERTER V«TH A 
NONZERO COi»MUTATION TIME 


3.1 INTRODUCTION 

This chapter presents the steady state analysis of the 
three phase squirrel cage induction motor fed by a current 
source inverter with a nonzero commutation time. It has bee 
shown in [l] that the rise or fall of the line currents of 
the motor during the commutation of the invertor can bo appi 
mated by a function proportional to I^^ cos m^t, where is 
parameter dependent on the leakage inductance of the motor a 
the value of the commutating capacitor of the invertor. It , 
also been shown in [l] that the time duration of the commuta' 
tion, denoted by t^ is given by 

t„ = Tt/Zw (3.1) 

c c 

The expressions for this assumed currents waveform have boon 
given in Table 3.1. Fig. 3.1 shows the waveshapes of those 
lino currents . 

The dq model of tho induction motor ^ with axis attached 
the stator (Sec. l.l ), is used for the analysis. The 
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transf orinod dq currents, using equations (4,39), for the 
three intervals of the inverter period aro also listed in 
Table 3.1 and ’I‘ in Table 3.1, is defined by equation (2, 
The Vi/avoforms of the dq currents aro given in the figure C 

The v/avoforra of the rotor current is obtained by the 
solution of equations (1 and the phase voltages of th 

motor aro obtained through equations (1,7). 

The analysis has boon done on the similar linos as th. 
followed in Chapter 2. The difference in this caso being ■ 
different expressions of the line currents. Sections 3.2 ■ 
3.5 correspond to the Sections 2,2 to 2.5 respectively. 
Section 3.2 deals with the time domain analysis, considorir 
both the stationary and the rotating rotor cases soparatel'v 
Section 3.3 deals with the frequency domain analysis and in 
Section 3.4 the voltage expressions at the stator terminals 
have boon obtained. In Section 3.5 the caso of the inducti 
motor introduced in Section 2.5 has boon considered. 

In this chapter, wo have an additional parameter j 

the value of increases, tho current waveshape tends to tl' 
ideal quasi square waveshape, corresponding to Fig, 2.1. 11 

has also been shown in this chapter that as tends towards 
infinity, the rotor currents and stator yoltago expressions 
this chapter approach those obtained in Chapter 2, 



3.2 MALYSIS OF INDUCTION MOTOR FOR ROTOR CURRENTS THROUGl 

TII^E DOIAAIN 

Tho analysis proceeds in similar lines as done in Sec, 
The equation of tho motor porformanco (1.6) is solved. The 
arc solved for three intervals (1,11,111) only as tho state 
waveform is symmetric about wt = it. Each interval is consi 
red in two portions. One corresponds to the commutation po 
and tho other to tho remaining interval period. So the fir 
interval corresponds to the time from start of interval, t^^ 
less then time of commutation, t^. This analysis has bee 
done in [12 ] . The results are reproduced here for all the 
three intervals. It should be noted that the symbols of in. 
tial currents correspond to equation (2.2), 

3.2.1 Case of stationary rotor 

The expressions of the rotor currents, as taken from [J 
for intervals I, II and III have been given below. 

3. 7. 1.1 Interval I ( o < wt Tt/S) 
for o < t < t 

N2 = ^do C3.2) 

iq 2 = ^^qo'^^^^1^ exp("at) ■-Y'3KjL[cos({o^t)~C2 sin(a)ct)] (3.3) 





3.4 

.1 

qo 

~f SICK, 

= [exp(-aTj » l/2) 

(3.5) 


- ivb>^l/(a^+a)?) L 22 

(3.6) 

C 2 

= a/w^ 

(3.7) 

K 2 

= [I+C 2 exp(at^)] exp(-aT^^) 

(3.8) 

T. 

= interval period = %/3ui 

(3.9) 


and * a’ , 'C^’ are as defined in equations (2.15) and (2.23) 
respectively, 

for t < t < T = ii/3w 

•*»-* »=*** 

^d2 ~ ^do (3.10) 

iq2 " ^q2 exp(~a( t-t^) ) (3.11) 

where 

’-q2 = (3.12) 

3. 2. 1.2 Interval II (-3 i i 3^ ° S . “*2 i 

/ 

for o < t < t 

msrs^ Vj» 

3Kj, 3 

^d2 = ^^do “ “T") exp(~at2)+ 5 K^[cos(a)^t2)~C2 sin(o>^t2)] 

(3.13) 



3.5 


f 3Kj; f 3K, 

^q2 " Gxp(~a1:2)- [cos(a>^t2) ~C2 sin(a)^t2)] 

(3.14) 


wnoro 


fSKpK 


■qo 


1 [1 . j 


3K2Kj_ 

^do [e^P(“aT^)] 


(3.15) 

(3.16) 


for < t2 < T^, i.e.((^ -i- t^) < wt < j£) 


^d2 = ^2 G^P('-a(t2 ~ t^)) 


(3.17) 


■q2 = ^q2 c^xp( “a( t 2 -t^) ) 


(3.18) 


v'hora 


1^2 = (I 


3K, 


do 


3 K 1 C 2 


— ) exp(-at^) “ — g- 


(3.19) 


fSK, f3K,Co 

^q 2 " ^^qo “ 2 *^^ Gxp(-at^) + 

3. 2. 1.3 Interval III (|2- < wt < TT i.e. 0 < 0)13 < |) 


(3.20) 


for o CO 13 < i.c. ^ cot j< (^^- + t^) 


■ 2 % 


3 3 

^d2 ” ^^do“ S' exp("'at3) + Kj_[ 003(00^13) -C2 sin(to^t3)] 

(3.21) 
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iq2 = (Iqo - ^ expC-atg) + ^ Kj_[cos(M^t 3 )-C 2 slnCoj^tj) 

(3.22 

where 


^do " I %7“ 

4 


Iqo ^ [expC-aip + l] 

(3.24 

for t_^ < mto < § 

C — 3 “■ 3 


^d 2 = ^^^2 ®xp(-a(t 3 -t^)) 

(3.25! 

iq 2 “ ^q 2 e^P("’a('t3“'tc) ) 

(3.26) 

'vVhero 


^d 2 = - 1 h) oxp(-atp - § Kj^Cg 

(3.27) 

^i 2 = (Iqo - '2 h’ oxpC-atp - ^ K 3 C 2 

(3.28) 


Tho expressions of the initial and final values of currents fo 
each intervals have been summarised in Table 3.2. It can bo 
scon from tho table that tho expressions for tho final values 
of currents in an interval are the same as the initial values 
of currents in next interval. Hence, tho rotor currents are 
continuous . 
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3.2,2 Case of rotating rotor 

The expressions for quasi rotor currents because of th« 
stator waveform of Fig. 3.1, have reproduced in this section 
from [ 12 ]. Actual rotor currents can be computed from these 
quasi rotor currents by using equations (2.44) and (2.45). Tl’ 
symbols used below correspond to equation (2.58). 

3. 2. 2.1 Interval I (0 < wt ;< ^) 
for o < t < t 

c 

X^(t) = exp(-“at) sin(m^t)-s-B^^ cos(w^t)] “• -f 

-}- sin(Wj,t)+K 2 cos(u)^t) ]-i-Kp^w^[c 0 ^ cosa)^t'-'K 2 sinw 

(3,29), 


whore 


Bol 

— V * 0 

- yvj; Q -r Uj 

(3.30) 


= 

(3.31) 

A , 
oi 


(3.32) 

2 

= (m^-fa^"~a)^)/ 2 a 

(3.33) 

K 3 

= i/(w'Ak|) 

1 

(3.34) 

>^4 

= ^ K^l (r3w^a)/(a^-to^) 

(3.35a' 
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K, 


'^3/^22 


(3.35b) 


K 5 = f 3K^IK3/2a 

Xio = (»rX20-2Xl0 + Kj^I 

and 


(3.36) 

(3.37) 


X2(t) sxp(»-at) [a ^2 sin(a)^t)-i*B^2 cos(a)^t)] + 


whers 


(3.38) 


'^^02 " ^0 ^3 
^3 “ 

^o2 = 


(3.39) 

(3.40) 

(3.41) 


Kg, = f3K3_I(a) -f3a)/2( a2+a)2) 

X" 


^20 "^r^lO ~ ” 5" ^1 


K,I 


^20 ^10 ®^P^Gssion is given by 


: X20 I ' i/2"(y4Y2-Yiy3) - +(Y 


'4"2 ^1^3" ~ 2 


I =^10 j pi 


S- -(Y4Y2-AY3)i 


(3.42) 


(3.43) 


aYa+YiY^), jz 


i z. 


(3.44) 
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where 


~ 

(3.45) 

Y, = 3Xp{-at ) sin a) t 
c re 

(3.46) 

Yj, = exp(-at ) cos w t 
c re 

(3.47) 

•/3 = oxp(.-a(T^-tp) sin[Mj.(T<,-tp] 

(3.48) 

Y 4 = exp[- 3 (T^-t )] cos[a) (T -t )] 

X c c 

(3.49) 


(3.50) 

^2 = Vl’^^3''^2'^^6"'^5^rS 

(3.51) 

C 4 = [ 033 - 1 K^I+K^w^to^l/w^ 

(3.52) 


(3.53) 

C 5 = [~ '^1 K^I+K4a-f-w^Z2]Ao^ 

(3.54) 

Z 4 = ^3V^4(Z2“K6)+K^ 

(3.55) 

C 5 = [“ZiW^~ 1 K 3 ^I-K^a)]/co^ 

(3.56) 
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The order of calculations of the constants for calculation of 

•^ 2.0 ^1 ^ 6 * ^1 ^ 4 -’ ^1 ^ 4 -* ^ 1 ^ ^ 5 ’ ^ 6 * 

and Z/,. 

for t^^ < t < I 

X^(t) = oxp[-a(t-t^)] sin(w^( t-^t^) + cos(w^( t-t^) 3 -K 

(3.57 

X2(t) = exp[~a( t-»t^) ] [0^2 sin(u 3 ^( t-t^^) )+Dq2 cos(w^( t-t^) )+K^ 

(3.58’, 


where 



°ol 

11 

X 

(3.59) 

^ol 

= [X2 + aD^j^ ]/Uj, 

(3.60) 


= Wl+^10^2«1 

(3.61) 

a 

= - aXi - K^I 

(3.62) 

^2 


(3.63) 

”^02 

11 

1 

X 

(3.64) 

^o2 

= [^ + 80^2 j/“r 

(3.65) 

^2 

“ -“r^i - =4 - 1 Kil 

(3.66) 
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3, 2.2. 2 Interval II (i.e. (i.e. o < o)t 2 < 

for o < to < 

~ 2 c 

Xj_(t) = exp(<-at2) [Aq 3 si5^(^^3-'t2)+Bj^3 cos(o3^t2)] + 

"-K 3 +Kq[w^ sin(ia^t 2 )+K 2 cos(w^t 2 )] + 

+ cos(w^t2)“-K2 sin(u)^t2)3 (3.67) 


X2(t) = exp(~at2) sin(o)^t2)+B^^ cos(w^t2)3 + 

"i" s in ( ( 0 ^ tp ) tK .2 cos(ci)^t 2^3 ■!" 

- K3^[ai^ cos(w^t2) - K 2 sin(to^t2)3 (3.68) 


where 





®o 2 

= >^10 ^7 " ^ 8 -^ “ 



(3.69) 

CO 

0 

< 

= [Xj_Q+aBQ 3 ~KQa)^+KgK2a)c3/w3, 



(3.70) 

^10 

= ( 0 ^X 20 - aX^o - ^ K^I 



(3.71) 

K? 

= (f3IK3a)/(a^+a)^) 



(3.72) 

^8 

= [(VSI K^a/2)-(3I Kjl u)j,/2)3 

(K 3 / 2 a) 


(3.73) 


= f3I K 3/4 


\ 

(3.74) 
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®o4 

^0“^'^10‘^^1 1^2 "^12“c 

(3.75) 

^4 

[X2o+aBo4+KiiU)^-K^2^2‘^c^ /'^r 

(3,76) 

^20 

3 

— '"^co^Xj^O ^^20 ** ^ 

(3.77) 

0 

t — i 

y: 

= ^31 K,m^/(a^+wJ) 

X X T 

(3.78) 

Kll 

= [(3Kj_ Ia/2)+(f3Kj_w^ l/2)](K3/2a) 

(3.79) 

=^12 

= 31 KgAa 

(3.80) 


The values of and X 2 Q is determined from matrix equation 


•20 

1 

1 

i I” 1 

, j ^ . 

[ 

(Y4Y2-Y;^Y3) 

—1 

- V3 H- (Y3Y2+Y4YJ,) 1 ■ 

00 

N 

1 

aoj 

j 

, f3 

i 2 “ 

(Y3Y2+Y4Y2) 

1 " ^’^4^2“'^ 1^3^ : 

( 

i 

j ^7 


(3.81) 


The constants ai'e to be solved in following order to 

Ky to ■*'^ 2 . 2 * ^1 A:* ^7 ^10^ ^ ^6^ ^11^ ^12^ ^7 

where 

KglC^ - (3.82) 

Cg = Caa^-(f3K^ l/2)-KQm^+KgK2m^]/(0^ 


(3.83) 
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^9 = (3.84) 

Cio= [-(-3Kj_ l/2)+aCg+K^jL^^-.Kj_2K2^c^/‘^r (3.85) 

^5 ^8^1'^^^2’’^7'^^8^c"^9^l2 (3.86) 

“6 ^10^l'^'^9^2^^10"^ll^c'^^12^2 (3.87) 

^11^ Kj_I-i-K^a]a)^ (3.88) 

'^12'^ (3.89) 

^7 ^11^3'^(25+^7^^4“'^7 (3.90) 

28 "" ‘^12^3'^^^6“^10^^V^'^10 (3.91) 

For t^ < cot 2 < I 

/l(t2) = exp(-a( t2-t^) ) [0^3 sin(w^(w^( t2-tp) ) + 

+ D_^3 cos(M^(t2-t^))]-K7 (3.92) 

^2^ ''-2^ “ expC-aCtj-tj,) ) [C^ sin “ 3 ,(t 2 "*c'^ ■*■ 

+ cos(Uj.(t2-tj,))] + Kj^g (3.93) 

where 

Do3 = 


(3.94) 
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n 

'^o3 

= [Xi + aD^g] 

(3.9b; 

^o4 

= X4 - K^O 

(3.96) 

Co4 

= [X^ + aDjj4 ]/Uj, 

(3.97) 


X”:, and can be obtained from equations (3.67) and (3,68) by 

• « 

substituting t^ = t^. X£ and can be obtained from equatio 
(2,50) and (2.51). 

3.2.2.^j Interval III (i.e, ^ tot ^ 7t or o ^ ^^3 S. 

for o < t _< t|^ 

X^( t) = exp(“at^) [Aq 5 sin(a>^t2)+B^^ cos(a)^t2)] + 

%3’‘%4^“c sin(to^t3)+K2 cos(a)^t2)] + 

“ slnioi^t^)] (3.98) 

X 2 (t) = exp(-at^) sin(a)^t 2 )+B^^ cos(a)^t 2 )] + 

**■ %6'‘'%7^^c sin(m^t3)-hK2 cos(w^t3)] + 

- fS %5[w^ cos(w^t3) 4 - K 2 sin(w^t3)] (3.99) 

where 

^o5 ^ ^10 ~ ^13 ^14 ^15 


(3.100) 
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"^05 (3,101) 

^10 ~ ^r^O ” ^^10 ~ (3.102) 

~ [ir3Kj|^l/2j (V*3{i)^“a) /( a^-Hi)^) (3.103) 

%4 "" [(f3K^I a/2)+(3K3_Ita^/2)] (K3/2a) (3.104) 

^15 “ f3K^Iw^K3/4a (3.105) 

®o6 " ^0~^l6"^17^^'^''^^^15^c (3.106) 

^o6 = ^^0^®®o6'“^^17‘^c“'^^-'^15^^2'^c^/''^r (3.107) 

« 

^20 = ^10 ■ (3.108) 

“ (V’3K2l/2)(iaj.+f3a)/(a^+u^) (3.109) 

Kj^7 = [(T3K^u>^ I/2)-(3K3^aI/2)](K3/2a) (3.111) 

For the values of ^10* have 


^20 


r* 

J “ ^^4'^2"'^l'^3^ - ^5 + ('^3'^2'*'^4^1^ 

-1 

j 

^12 

>^10 


- (YjYj+Y^Y^) § - (Y^Yg-y^Yj) 

( 

1 

I 

N 

i 


( 3 . 111 ) 
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i he constants are to be evaluated in the following order 


Ki to K3, I<33 to K37, Y3 to Y4, C33 to 
and Z-,^, where 




n 

JL 

= ‘"K^3+K^4‘<2+^15^c 


(3.112) 

^14 

= [ “'r3Kj_I+aCj_3+K^4W^--K^^K2W^ ] /w^ 


(3,113) 

^15 



(3.114) 

^16 

= [aCj_5-K^^(^2^f3Ki5K2U)c]/w^ 

- 

(3.115) 

“9 

^14^1'^^1 13“%4‘^ c'^^1 5^ 


(3.116) 

^10 

^16'^l'^^15^2'^^l6^^17^c'^^^^'^15^2 


(3.117) 

^17 

~ {- Y"3KiI/2 ) — Kj^sO 3 /Wj, 


(3.118) 

^13 

~ L”Z9^3,+ (3K^l/2)-K^^on3/w^ 


(3.119) 

^11 

= C;^7V^4^Z9-Ki3)-fK3_3 


(3.120) 

^12 

= Ci3^3+(2:io~%6)VKi6 


(3.121) 

For 

■^c 1 "^3 1 1 




X^(t) = exp(~a(t3^.t^))[C^^ "^r^ ^ ^-^^13 


(3.122) 
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X^Ct) = exp(-a{t3«»t^)[C^^ sin(w^( 13-!^) )-5-D^^ cos w^( 13-!^) )J+K 

(3,123 


where 

“os = - Kis 

'^' o 5 ~ ^^1 ^ 

^06 ^ “ '^16 

^06 = ^ 


(3.124 

(3.125 

(3.126 

(3.127 


The expressions for Xj^ and X^ can be obtained from equations 

* 

(3.98) and (3.99) respectively by substituting ^3 = 1^, X]_ and 
X2 can be obtained from equations (2,50) and (2.51) by sub- 
stituting for X|^ and X^, A computer program has been develope' 
to calculate initial and final values of currents for those 
intervals in the rotating rotor case and the boundary condi- 
tions are found to be matching (Sec. 3,5). 

3.2.3 Case of to — v <» 

As 00^ increases, the stator current vraveform of this 
chapter, (Fig. 3,1), tends towards the v/aveform considered 
in Chapter 2, (Fig. 2.1). In this section it has been shown 
that as to^ tends towards infinity, the expressions for rotor 
currents obtained in Sec. 3.2,1 and 3.2,2 reduce to that 
obtained in Sees. 2.2.1 and 2.2.2. This reduction is shown 



3.18 


for interval I for both tho sections. The results for other 
intervals can bo shown by proceeding on similar lines, 

3. 2. 3.1 Stationary rotor case 

As u) oo wo note from equation (3.1) 

t^ = 0 (3.128) 

So tho interval o < t < t of Sec. 3. 2, 1.1 reduce to a point 
t = 0 and the period t < t < to tho remaining period of 
interval I except t = Oy as is tho case in Sec. 2, 2. 1.1, 

Also, as ' oo 
c 


K = JiiL- = K» [From (3.6), (2.13)] (3.129) 

j. ^22 

C 2 = 0 [From (3.7)] (3.130) 


K 2 = oxp(*'aT^) = [From (3.3) , (2.23) , (3.130)] 

(3.131) 


Thus from (3.129) and (3.131) in equation (3.4) reduce to 

that in equation (2.21) and of equation (3.5) to that in 

qo 

equation (2,22). Also 1^2 equation (3,12) reduce to 
Iq2 ~ ^^qo snd this equation is similar to (2.12). 

The expressions of currents in Sec, 3, 2. 1,1 become 


^d2 = ^do 

S 2 “ ^q2 


[From cqn. (3.10)] 


(3.132) 


[From eqn, (3.11)] 


(3.133) 



3.19 


ThosG expressions are equivalent to equations (2.16) and 
(2.17) os (3,133) is already shown to reduce to that in 

(2,12) and 1^2 oqn. (2.10) in equation (2.16). Thus 

the expressions obtained in 3, 2, 1,1 reduce to that obtained in 

2. 2. 1.1 as approaches infinity, 

3. 2. 3. 2 Rotating rotor case 

Hero it will be shown that for co approaching infinity, 
the rotor currents expressions of section 3.2,2, 1 reduce to 
that obtained in Section 2, 2, 2,1. It has already boon shovm 
in 3. 2. 3.1 that sub regions of the interval taken in 3.2.2. 1 
reduce to that taken in 2.2 .2.1 as --- «>, 

For t = 0 

So, eqn. (3,29) is for a point t = o. Thus from (3.29) 
X^(t=o) - ” ^^4 + ^^^2 (3.134) 

Substituting for from (3.30) in (3.134) 

X^(t=o) = X^Q+C^-K^'fK 5 K 2 a + (3.135) 

Substituting for from (3,31) in (3.135) 

X^(t=o). = (3,136) 

as should bo 

Similarly oqn. (3,38) reduce to 

^2(t~o) = • 
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For and X 2 Q 


As «>, since t = 0 , we have 

= O [from (3.46)] 

= 1 [from (3.47)] 

Y^ = exp(-aT^) sin [from (3.48)] 

Y 4 = exp(”aT^) cos [from (3.49)] 

Cf = [from (3.31), (3.33), (3.34), (3.36)] 

^ 1=0 [from (3.45) , (3.36) ,(3.34)] 

C 2 = [- KiI+aK 4 ]/Wr [from (3.50)] 

= ~K^ [from (3.40)] 

C 4 = [-*K^a - K^l]/w^ [from (3.52)] 

Z 2 = 0 [from (3.51) ,(3.136) to (3,144)] 

^5 ~ Kj_I+K 4 a]/a)^ [from (3,55)] 

= [- § K^I ~ K^a]/w^ [from (3,56)] 

Z 3 = C^Y^ + K^Y^ • [from (3.53) , (3.141)] 

^4 “ ^3'"6“'^4^6'^^% [from (3.56) , (3.145) ] 


(3.136) 

(3.137) 

(3.138) 

(3.139) 

(3.140) 

(3.141) 

(3.142) 

(3.143) 

(3.144) 

(3.145) 

(3.146) 

(3.147) 

(3.143) 

(3.149) 
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and from (3.44) using (3.136) and (3.137) 

(3.150) 

Thus it is observed that and of eqn . (3.146) and eqn. 
( 3 . 147 ) correspond to and obtained in equations (2.72) 
and (2.75). Also Y 3 and of Sec. 3. 2, 2.1 in reduced form 
for w in equations (3.138) and (3.139) correspond to 

and Y^ of Section 2. 2. 2.1, i.e., equations (2.68) and (2.69). 
Thus Z 3 and 2^ from (3.148) and (3.149) are seen equivalent to 
and Z}.^ respectively, as in equations (2.71) and (2.74). 

Thus equation (3,150) give the same value of X 2 Q and X^q 
give the same value of X 2 Q and 21 s obtained from equation 

(2.76), implying the similar initial conditions. 

For t > o 

The equations (3,57) and (3.58), using (3.123) can be 
seen to reduce as 

X.,(t) == exp(-at) sin “r* " °ol COS ^ ^ 3 ' C ^ 

X^Ct) := exp(*-at) [C^2 ^o2 (3.152) 

It can also be shown from equations (3.59) to (3.66) that as 

(jQ - ' CO 

c 
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= ><10 

(3.153) 

'/ t — V 

“2 ^20 

(3.154) 

><1 “ “rX20”>*Xl0 ■ 

(3.155) 

><2 = ^rXl0-^X20 - 2 

(3.156) 

^10 ^4 

(3.157) 

=01= + =‘Dod/“x 

(3.158) 

=02= ><20 “ ‘<6 

(3.159) 

=02= ^= 02 ] /“r 

(3.160) 


It is clear from equations (2,59), (2,60), (2.61), (2,66) and 
the above equations that 

0^2= a„2 (3.161) 


and 

=i2= =02 


(3.162) 


Thus equations (3.151) and (3,152) correspond to equations 
(2.56) and (2.62). Thus, as w -- «>, the results obtained in 

C 

Sac, 3, 2. 2.1 reduce to that obtained in Sec, 2. 2. 2,1, as it 
should be. 
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3.3 ANALYSIS OF INDUCTION MOTOR FOR ROTOR CURRENTS THROUGH 
FREQUENCY DOMAIN 


ThG analysis to obtain the rotor currents through froquenc' 
douains has two steps, as is done in Soc, 2.3. First is the 
calculation of the harmonic components of the stator current. 
Then the harmonics of quasi rotor current, and X2 are com- 
puted by the use of equations (2.52) and (2.53), Tho harmonics 
of actual rotor currents can be obtained from equations (2.42) 
and (2.43) after wo have computed the harmonics of stator and 
quasi rotor currents. 



For the harmonic components of the stator current, 
and i^2^ can be written in Fourier series as 

00 OD 


Fig .3.1, 




cos nmt + 


E 

n=l 


d^ sin ntot 


(3.163) 


i 1 = I C cos nwt + 
^ n=l q 


00 


E 

n=l 



sin ntot 


(3.164) 


From Fig, 3.1, it can bo said that for even values of n, the 

constants C„ , C„ , d„ and d are zero because of mirror image 
^d "^q 

symmetry. For odd, values of n. 



[-sin 



2 _ 

no) 


1 

(nw+m^7 


1 

Inw-a)^) 


+ 


. 1 - 1 





)+cos ( 


nut) 



+ |2); ] 

(3.165) 
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°n^ - T ^ TSS^J ■ •nsmop'. == *"'"^23; 


nTO\ . 


i. < Sn-n; '. ' n 


r- - 2 nz . 

^n “ “fa 3 ^n, 

q » c 




r.Tt 


'n T3 '-n. 


(3.166) 

(3.167) 

(3.168) 


It can be seen from equations (3.166) to (3.168) that for n to 
be multiple of 3, these constants are zero. Thus from noWj. 
when we refer to hanBonics v/e mean onl y odd and non tr ipl ents. 

Thus t writing the mth harmonic component of i^^. iql 
as 


^dlm = Im, cos(nMt + p ) 

U G 

^qlm ~ ^raq + ^mq^ 


We note from equations (3.165) to (3.168) that 


1 

= jl’ 

[ = 

U’ 1 











= 

(Pm - 
. 

rtx 

j) 

for 

m 

= 1,7,13, 

19 ... 



ITv 

for 

ra 

= 5,11,17 

, 23 , ... 

q 

d 





If we 

define 

•p‘ 

as 

in 

equation 

(2.123) then 


(3.169) 

(3.170) 


(3.171) 

(3.172) 


(3.173) 
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~c.lm m 


= cos(imot + 


(3.174) 


qlff. = -m 5 ) 


(3.175) 


where 




2 ^1/2 


(3.176) 




(3.177) 


It should be remembered here, that m refers only to odd and non 


:iplent value. Solving (2,52) and (2.53) for mth harmonic of 


X-, and X 2 respectively, as done in Sec. 2.3 we can again visua- 


lise the picture as in Fig. 2.3. For this case we have 


K^a(»l)P 


I’ sin(rnwt - a. ) + 


+ (m -(-1)^ mm) cos(in£ot + p 

Af - m^ 


(3.178) 


K,al 


1® m 


sin(mwt + - ctj) + 


K T' ■ 

+ — (m =(-1)^ mm) sin(mmt4-6 -«_,) 

J(\rn JL i 


(3.179) 


where Aj and are given by equations (2.129) and (2.130). It 
can bo noted from equations (3.170) and (3.179) that the 


magnitude of harmonics and X 2 jjj are equal and ihey have a 
phase difference of ti/2. In Sec. 3.5, harmonics of X^ and X 2 


have boon computed for various values of m and m^ for a case oi 

W * 


induction motor. 
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CaSG of w 

As Ci>^ incroasos, the stator current waveform of this 
chapxar (Fig. 3.1) toads towards the v/avoform considorod in 
Cheptor 2, (Fig, 2.1), So, in the limit of to tending towards 

V# 

infinity, the result of this section tend to that obtained 
in See. 2.3, as is shown bc-lov; 



(3.180) 

(3.131) 


X • nj • 

Lt C = (3,132) 

m - c= i^d “^d 
c 

Thus, from equations (2.126), (3,176), (3.180) and (3.182) 


Lt 





From equations (3.177) and (3.180) 


(3.183) 


.6,, = 0 (3.184) 

i^d 

Substituting for I* and 0-, from equations (3.183) and (3,184) 

ra 

in (3,178) and (3.179) we see that those expressions of har- 
monics of quasi rotor currents reduce to that obtained in 
3oc. 2,3, equations (2.127) and (2.131). 
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3.4 VOLTAGE WAVEFORM AT THE STATOR TERMINALS 


The expressions for the voltage waveform at the stator 
terminals for the stator dq currents as in (Fig. 3.2), can be 
obtained by proceeding similar to Sec. 2.4. For this analysis 
each interval is divided into three regions. They correspond 
to time (t- < t-), (t. = t_) and the remaining period of the 
interval. In this section the cases of the stationary and the 
rotating rotor have been dealt with separately. 

3.4,1 For the stationary rotor 


The equations for stators voltages for this case of 
stationary rotor have been computed in Sec. 2.4.1. They are 
equations (2.137) and (2.133), i.e.. 


^dl ^l^dl *“ ^d2 ^^1 *■ 

m2 

Vql ^ ^l^ql ^ ^q2 "*■ ^^11 " I^^P^ql 


(3.185) 

(3.186) 


Since i^j^ and i^^^ are continuous, (Fig, 3.2), there is no 
impulse in this case. 


At t^ = 0 



The jump in V^^ and V^^ at t^^ =0, given by 
can be computed from (3.185) and (3,186) 


A'V^l and 
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'dl 


= (L 


11 








‘dl' 


t.=o 



(3.187) 


(3.188) 


Table 3.1 can be t^sed to compute the changes in slopes at 
= o. From Tablo 3.1 it can be shown 

0 ^ 3 . 189 ) 

j 1 


'-\(pi 



t.=o 


0 


(3.190) 


where 1 = 1,^,3,, ... 


and so '^dl and continuous at t^ =0. 

For 0 < t. < t 
-1 c 

The equations (3.i85) and (3.186) can be used to compute 
the ex^sressions of the yoltages. Table 3.4 summarises the 
references to the expre^sioios of ^q2 

required for the computation of voltages. 

At t. = t 
1 . c 

• -^is=7Tir r «. t 

It is seen from (Fig. 3.2) that there is sudden chance of 
dopes of i^j^ and i^-, at this point. Hence ^ there is a jump 
in the voltagxo expressions at t* = t^, as can bo seen from 
equations (3.137) and (3.188). 
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Frora TeJ^lo 3.1, using = ofn" f^rom equation (3.1), wo have 

V,# ^ 

c 


dl^ j 

(t=tp “ ° 

(3.191) 

al> 

(t=tp =^31^0 

(3.192) 

dl)| 

kt2=t^)“ " 1“ “c 

(3.193) 

) i 
ql^l 

!(t„=t )°'' % 

1 ^ ^ 

(3.194) 

di^ 

1 

_ 31 

(t3=t^) 2 

(3.195) 

■ql^ 

i 

1 f3I ,, 

|(t3=tc)=‘^ 2“"^c 

(3.196) 


Equations (3.i9l)to (3.196) can bo used in equations (3.137) 
and (3.138) to cOiTiputo the jumps in and at ~ "^c* 


For < t. < 

c 1 c 


During this interval i^j^ and i^^j^ are constant. Hence, 
equations (3.185) and (3,136) roduco to 


''di = id2 




(3.197) 

(3.198) 


Tabic 3.3 gives tlio roforonces for the computation of the 
voltages . 
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3.4.2 For the rotating rotor 


From equations (2.148) and (2.149), using equations 
(2.42) p (2.43) to elliminate for ±^2 ^q2 have 


V^l — Mayv^ ~ (3.199) 

'■'ql = ■‘^^1 - + '*>A + V^ql (3-200) 


where and are defined in equations (2.140) and (2.143) 
respectively. 


Since i^^ and i^^^ are continuous. Fig. 3.2, 
are also continuous, equation (2.45) and (2.46), 
impulse in the expressions of and Vq 3 _» 


and 

there are no 


At - 0 


It has been seen in Section 3.4.1 that at t^ = 0, the 


changes in and ’pi^^^’ are zero and hence and V 


ql 


is continuous at t. = o. 


For o < t- < t 

1 c 


The expressions for and X^ obtained in Sec. 3.2,2 to- 
gether with Table 3.1 is used in equations (3.199) and (3.200) 
to obtain the expressions for voltages. Table 3,5 lists the 
roferencc-s to equation numbers for Xj^ and X 2 - 
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Fo r tj_ = 

From equations (2.199) and (2.200) because 
X2 are continuous 


'^dl^^i~^c^ “ ■^1 ''P^m ) 






t.=t 
1 c 


%=^c 


(3.201) 

(3.202) 


The e quations (3.191) to (3.196) can be used in equations 
(3.201) and 9(3.202) to obtain these discontinuity of voltage 


For t < t. < 

c 1 c 


In this region i^^j^ and i^j^ are constant and hence equa* 
tions (2.199) and (2.200) reduce to 


^^dl ^ ^ ^dl ** (3.203) 

\-ll ^ ■^‘2 ^ql " ^"'''^^^1 ^ (3.204) 

Table 3.4 lists the references to the equations which give th« 
expressions for ^i>^2’^dl ^ql their region, and 

can be used to compute v^j^ and 

3,5 CALCULATIONS FOR A KNOFJN fvlACHIjNIE PARAMETERS 

Using the results of Sec. 3.2 and Sec. 3.3 two computer 
programs are developed to compute rotor currents. The listinc 
are given in Appendix B. 



3.32 


Using those programs the data for Tables 3,5 and 3,6 arer 

obtained. It is seen from Table 3.5 that the I and I „ 

■ do qo 

values need to computed only for one interval. Table 3*6 
shows that the solutions obtained through both the programs 
using Sec, 3.3 and Sec* 3.4 results are similar. It is also 
seen from Table 3*6 that as w increases, the results approach 
to that obtained in Table 2,6 for ideal inverter current. 




Three- phase and two phase stator currents 
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of rotating rotor case 
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TABLE 3,5 


TiiTEHVAIj 

HU.viBE!? 


I 

II 


** P’UR XI** 
iMiTXAb FINAL 

value SfALUE 


FOR 

INITIAL 

VALUE 


X2** 

final 

i^AbUE 


FOR ir>C=b,A'4fc»;l«V. FKEO.=20,HZ|PPM*§00,?.^Ca500, 
i 3.2625U 0,69P4 -1,09295 -3,37006 

II C.6fe468 -2.57277 -3,37165 -2,28387 

III -2.57730 -3.26243 -2,27890 1.09291 


II 0,92027 -0,65605 

III -0,85950 -1,77973 

FOR ir>C*I.Af4FyINy. FREQ, = 15 


1.07492 0.41953 

«. 41779 -0.65515 

III -0,65713 -1.07490 

FOR I0C»6,AMPnNV. FREQ, =20 


,H7,?FPM=600,?^C=800, 
0.035P -1,52152 

-1,52376 -1.56027 

-1.55888 -0.03612 

,HZ?BPMs200.1«C=750. 

-0,13816 -0,99872 

-0.99998 -0,66266 

-0,86181 0,13816 

.HZ?PPM=400.?i<JC=700, 


III 


1.94711 

0,74558 

-1.20152 


0.74907 

-1.19805 

-1.94707 


-0,26321 

-1,817§3 

-1.55462 


- 1.81507 

-1,55661 

5.26320 



table 3,6 

♦ ♦ f ^ ijt j|! * if; ^ 



FRRQ, 330 , Hg St WCfi 800 
, 1,54655 85.64724 

- r:0i98| 22.87026 

. 0.01107 -4.92045 

0,00289 ,,,, -62, 90466 
0,00182 266,05194 

0.00056 193,78426 

0*000,35 «.g,.l48;236l8 : 

1.06010 99,05103 , 

■■^ 0,02200 ■• 64 ; 4 | 3 y 3 .:'’'‘ 

0 101285 S3, iff 85': • ^ 

-0:004.||,^.', .31-06987 , 


85.64456 
22,83890 
- 4,90670 
: - 62,91205 
266,12543 
, 194.26089 
, i 48 ,. lS 40 S 

DC s^ 4 .AWPiINV 
99 . 04 a |5 




—————FOP XI————— 

i/lA TIMS !)Q.’4ft|i^ VIA FREQ, DOMAIN 

'•lAdUTUofi phase MAGNllUOB PHASE 


FOR HPH s 100 ,?IUC = 6 ,AMP|Il»V. 


i 0 , 8047 b 

5 0.01918 

7 0 .<,= ()740 

il 0,00105 

U 0.00068 

17 0,00038 

19 0.00025 

FOR RPH = 500 .; I DC 
I 3.41081 

6 0 . 01 701 

7 0,00819 

11 0.00098 

13 0.00072 

17 0.00038 

19 0.00025 


7 , 3.12485 
- 31.11028 
- 84,51469 
134.53249 
42,42627 
261.08836 
191.31723 

,s , 6 .AHP 1 INV, 
98,77277 
- 31.27879 
- 84,35751 
134,71308 
^ 42,24791 ■. 
261,73957 
191.68243 


FREQ, a20,HZ & WC= 300, 

0 . 80476 73.32652 

1 . jmi - — 


.,60740 
0,00106 
0.00066 
0.00038 
0,00025 


- 31.10732 

- 84.54451 

134.53023 

42.55904 

260.81009 

191,13576 


FREQ, = 20 , HZ 6 »«C= 300 


OlOSl 

0,01698 

0.00820 

0,00100 

0.00071 

0,00037 

0,00026 


98 , 7/671 
- 31,24336 
- 64,44720 
134,49919 
^ 42,58531 
260.79729 
191.14775 
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■fAHfiF: 3,0 


FOR KP' 

1 
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U 

n 

17 

19 

FUft RP'J 
1 

6 

i 

11 

13 

17 

19 

FOR RPM 
i 
5 
7 

11 

13 

17 

19 

FOR RPM : 
1 
5 
7 

11 

13 

17 

19 

FOR RPM 

5 

7 

11 

13 

17 

19 

FOB RPII 
1 

' 5 

7 

11 
„ 19 


•'"-"*FuR XI-.-.—,-. 

‘fi i 1 1 i ilUii i-nhiit FAGi'JTTHiJF PhA^F 


™ lOt., ;?:oG 
< f , R \) 4 7 o 
u,(-j d 
0.' '>7 40 

0 « 0 1. 'j O I’i 

O.OOt Ju 


= ^6.A0P;Ii’«V, 
73.3748b 
-31.110?b 
“«4,bl-ib9 
131.53249 
. 42.42627 
261.0RB36 
191,31723 


F 


1.4 f ';R1 
0 , I 7 (; i 

0 .ooi;9d 

9 .Out 36 
9.00025 

= SOO.JtfJC 
1 • 3 io5H 

6.6 I 9B3 

6 . V' 1 lOB 

O. 60289 
6.6U1 82 

P. OOoSu 
0,06635 

= 2oo.;ir)c 

1 .06«,*it 2 
0.02199 
0.01286 
0.00 447 
0,06335 
0,60170 
0,00137 

= 560,, 'IOC ' 
3.45 353 
9,02313 
0,0158? 
0.00524 

0,00433 

0.oi)224 

0.00199 


= 500 
. 3.45,.,^ 
0,02343 
0.0158B 
0,00524 
0,00434 
0.00225 
0.00199 


IDC 


=^b.AMP;INV. 

98,77277 
-31,27879 
-R4, 35/51 
1 34.71308 
_42, 24791 
261,73957 
191,68243 

= 8.AMP;lWv, 
85.64456 
22.8 3890 
'-4,90470 
-62.91205 
266,12543 
194.26089 
148.15403 

= ^4,AMPjrIf»V, 
99,64895 
64,46595 
53.38764 
31.00 1 3B 
20.64493 
—2,42762 
-13,52101 

= «*A0P|1NV, 
122,/l66u 
90.60155 
90, tBOOl 
89,71415 
89,54617 
89,09412 ■ 
86,94177 

122178717^^' 
91.04879 
90,97852 
90,50157 
90.47295 . 
90,|li04 
90,29401 


0,804 76 
0,01917 
0 . 067*0 
• 0.00106 
0,0u068 
0,00038 
0,00025 

FRFQ. -20 
3.41051 
0.01698 
0.06820 
0.00100 
0.0607 I 
0.00037 
0,00026 


~20,il7 S yJC 


30 U 


73,32652 
-31,10732 
-84,5445} 
134,53(:.2 3 
42, 559 1' 4 
?60.Rlf;b« 

1 9l . 1 3576 

HZ 6 wC=: 300. 
98.7/671 
— 3 1 , 2 4 3 3 6 
-84.44720 
134,49919 
,‘l2,5a5il 
^00.79729 


191.14775 


FRF6, =36, HZ i, WC': 


' ^ 

1,34655 
0.01983 
0.01107 
0,00289 
0,06182 
0.00056 
0.00035 

FRFy, =15 
l.OoOlO 
0.02266 
0,01285 
0,00448 
0.003'35 
0.00170 
0.00137 

FHFO, =20. 
3 . 45320 ’ 
0.02347 
0,01586 
0.00526 
0,00433 
0.00226 
0,00198 


860 


85.64724 
22,87026 
-4,92045 
-o2 , 964 66 
2o6,05t 94 
193,78426 
148,23648 

,HZ 0 wc= lOon, 
99,65103 
64.48393 
53,38785 
31.0 u987 

-21^9^2 

-13,53293 

HZ 6 wC=99999. 
122.72068 
90,69405 
96.49350 
89.72814 
89,50972 
89.12122 
88,97076 


^8|Q. =20,H2 it «iC5:9y9999« 
3*|52^2 122,79196 


foim 

olooill 

0,00433 

^•§§226 

0.00198 


91,05105 

90.99247 

90.51266 

90,49623 

90,33283 

90.32490 
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computatio:-: 


OF Tris ELHCTRO:ilGKETIC TORQUE OF 
CURREilT FED INDUCTION MOTOR 


THE 


. 1 INTRODUCTION 


ThG expression for the instantaneous electromagnetic 
torquo in terms of dq stator and rotor currents is given by 
the equation (1.8). This equation, together with che expre- 
ssions of dq stator and rotor currents, has been used to 
compute the torque produced by the squirrel cage induction 
motor fed by a three phase current source inverter. 

In the previous chanters, both the time domain and the 
frequency domain expressions for the dq stator and rotor 
currents vvore obtained. The time domain expressions obtained 
for the dq stator and rotor currents can bo directly used in 
the equation (1.8) to obtain the instantaneous torque v/aveforin 

The torque spectra -can also be obtained using the harnioni; 
component values of the stator and rotor currents in equation . 
(1.8). The instantaneous torque expression is given by the 
forms, which are product of the stator and rotor instantaneous 
currents. Thus if a nth harmonic of the stator current and an: 
rich haraionic of the rotor current, referred to the stator is ; 
assumod then these will produce (flH-n)th and the (m“n)th harmoni 



4,2 


component of tho torque. It has been shown in this chapter 
that the contribution to tha torque is only because of one 
of these- terms. Further it has been shown that the torque 
harmonic components are multiples of three for the case of ge 
X'al stator current. For the case of motor fed by the curren 
source inverter, tha torque harmonic frequencies are multiple 
of six times the inverter frequency. 

A torque harmonic has contributions from the various 
combinations of the stator and rotor current harmonics. But 
as the frequency of the harmonic current increases the contri- 
bution of these to the torque harmonic goes on decreasing. Ii 
face the contribution due to two dominant harmonics of curren' 
is seen to aporoximate fairly v;ell to the actual value of the 
torque harmonic. 

Sect ion 4.2 the nature of the harmonic components of 
the torque in any arbitary reference frame is discussed. The 
torque harmonic frequencies are found to multiple of three 
times the fundamental frequency for the case of general statoi 
current. The general expression for the torque produced due 
to nth harmonic of stator current and m'th harmonic of rotor 
current, -vhen referred to the rotor, has been computed in this 
section. 

Soc tio n 4,3 gives the procedures for the computation of 
the torque. One uses tho instantaneous stator and rotor 
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current expressions to obtain the electromagnetic torque 
produced. In the other the harmonic components of the torque 
have been computed using the harmonic component values of tho 
stator and rotor current. 

In Sec bion 4 ,A the stator current v/aveform of Fig. 2.1 he 
oeen used to compute the torcfue produced, using both the pro- 
cedures of Section 4,3, The harmonic components of the instai 
tnneous torque v^/aveform obtained are calculated. These are 
compared vdth the torque harmonic values obtained directly 
through the other procedure. This verifies that both are 
procedures compute the same result. 

4.2 NATURE OF THE HAPJ IONIC COiMPOMENTS INI THE TORQUE 

In this section, the nature of the harmonic components 
in the torque produced by tho squirrel cage induction motor 
has been established. To start with a arbitary stator curreni: 
waveform is studied. The expression for torque harmonic pro- 
duced due to the interaction of nth harmonic of stator current 
and mth harmonic of the rotor current, when referred to the 
stator, is coiaputod, in a arbitary reference frame. 

It turns out that torque harmonic is independent of the 
ch.oico of the reference frame and the torque hari.ionics have 
frequencies which are multiples of three times the stator 
current frequency. 
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stator current of the motor '-vhen fed by the current 
source inverter is syrnetric about oat = ix. Also it does not 
have triplens. The nature of the toroue for such a case of 
the stator current is also discussed in this section. Torque 
harraonic freauencies are raulripies of six tiiiies the inverter 
frequency in this case. 

4.2.1 Torque Marmonic in an arbitary reference frame 


Tiro torque harmonic produced by the squirrel cage induc- 
tion motor fed by a three phase periodic nonsinusoidal current 
of frequency oa, has. been computed in this section. The 
stator current harmonics 'vill have the frequencies of the 
typo 'nca', whore n is a positive integer. The rotor currents 
referred to stator are computed after applying a transfer 
function (Fig. 2.2) to the stator currents. Thus the fre- 
quencies of tho rotor currents referred to' the stator are alst 
of tho type 'mw'. 


The torque produced due to interaction of a stator 
current of frequency ’no’ and a rotor current of frequency 
’noj’j vdnen referred to tho stator has boon studied in this 
seebion. The rotor is assumed to be rotating v;i th a constant 
speed in tho direction of field produced by fundamental 
stator current. Then tho harraonics v'ith value of * m’ of the 
type (3p-'l) produce the fluxj revolving in opposite direction 
to the rotor speed. Hero, p = 1,2, 3, 4 ... 



4.5 


Thus the actual rotor currents will have the frequencies of 
the type (noo + lo^) where positive sign is taken for the case 
v/hero ’m’ is of the type (3p"l), and negative otherwise. 

The expressions for the stator and rotor currents with a 
arbitary values of magnitudes and phases, I , a and a 





5 

X £ 

respectively, can b 

>e written 

as 



^al ~ ^s 

cos[nwt - 

a ] 


(4.1) 

^bl “ ^s 

cos[nCti>t • 

3 ^ 


(4.2) 

hi = h 

cosLn(cot • 



(4.3) 

ha = h 

cos[ma)t ;h 



(4.4) 

^b2 “ ^r 

ros[m(a>t • 

“ i^) ± 


(d.5) 

^c2 "" ^r 

cos[ra(tot • 



(4,6) 


where lo'wer sign is for the case when m is of the type 
(3p-l). Upper sign is fox' all other cases. 

The currents in an arbitary reference frame [3] are related to 
lino currents etc., from equations (1.1) to 

(l«4), as 
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icu = (2/3) 
ic,i = (2/3) 
ie2 = (2/3) 
i,,2 = (2/3) 


[iai cosQ cos(0 ~ |>^3) + cos(© + |2.) ] (^.7 

sin9 - sinCO—^’^) - sin(0 + 3 “)] (4.8 

[ia 2 cos|3 + cos(p - | 2 .) + cos(p -h | 2 l)] ( 4.9 

[“ia 2 " ho2 “ §^^"^c 2 3 ^')] ("I-IO 


P = 0 - 0 (4.11 

r 

0 ^ = C0^,t (4.12 


In above is the constant rotor speed. (d0/dt) is the speed 
ox the arbitary reference frame. 

The torque is related to dq currents from equation (1,3) as 


T — 


1 

o'- ql 


•d2 


■dl 


^ 2 ) 


{4.14; 


\/he::G 


I'v = iM(m/2)(P/2) , (4.15) 

-loro m is the number of phase and P the number of poles. 

Substituting for dq currents in equation (4.1^1) froia 
equation (4.7) to oqn. (4,10) and sinplifying gives 
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h 

= (‘''/9) U 


t) 

sin 

(^» 

Q) 

+ e, 

2(t) 

s' 

Ln(;3- 

.0„ |iE) 

+ 






+ e3( 

t) 

sin 

(p 

- 0 + 

2.% 

■ 3 

)] 





(4, 

,16) 

o/hs 

‘XQ 

















ei 

X 

(t) = 

(i 

al 

^a2 

T“ 

^bl 

^b2 

4* 

"Cl 

^c24 


(4, 

,17) 



^2 

(t) = 

(i 

al 

~b2 

*r 

^bl 

^2 

4- 

^Cl 



(4, 

,18) 



©Q 

(t) = 

(i 

a2 

^bl 

4- 

^b2 

^cl 

4- 

^c2 



(4, 

.19) 

Suh 

s titu ci 

.ng 

for i 

al^ 

t) , 

^a: 

0 ( t 

) ia equal 

bion 

(4.17) 

f rofii 



equ 

,? tions 

(^. 

1) to 

(4 

.6) 

» Q 

ive 

s 









e^^vt) = igl^[ cos( ncot-Kg) cos(mt.ot + - a^) + 

+ cos(n(cot“ 'i^)-=-a ) cos(iTi(wt“ ■^') + w t - a ) + 

-i- cos(n(wt" cos(iii(cot“ + oJ^t-cc^)] 

This implies 

e^(t) = (lgI^/2)[cos [ (n+m)ajt“(ag+a^)+a)^t] + 

O-TT yi 

+ cos [ (n+m)(cot“- ^’)+ w^t~(ag+a^)] +cos[ (n+m) (wt- 3™) 

± [(n-rp.)a)t + co^t-( a^-a^) ] + 

2.71 *=- 

+ cos r (n-in) (cot- v“) - u)_t-( rc^-a^) 1 + 

+ cos [ (n-m) (cot- ^ 


( 4 . 20 ) 



4.8 


ihis can be shov/n to siiiiplify to 


= (3I,I.v''z) [pg cos [ (n4-n)(jt+w_t“(a^-!-c.^) ] + 
p,, cos [ (n“:a)wt T - (cc^-ctj ] ] 


(4 .21) 


P< 


if (n-rra) is of type ’3p’ 


0 


= 1 if (n-m) is of type ’3p’ 


>therv;ise 


(4.22) 


(4.23) 


Here p is any integer. 

Sirdlarly fro'-n equations (4.18) and (4.19) s it can bo shown 


C ilS 


&o('t) = (31 I^,/2) [p cos [(n+m)cot - + to t-«( cc t-c ) ] + 


'2 


s r' 


+ cos [(n“m)o)t + + co.^t - (a^-o:^)] ] (4.24) 

,( t) - (3IgI^/2) [pg cos [ (n-}-m)u)t+ 


2%m t , , •- f ^ , r., \1 , 

" L,“-( a^+o;^) J + 


1- p., cos [(n-r.)cot - + to.^t--(c^-G )] ] (4.25) 

o J- X* 


The expressions of Oi^'^t), ^tcom equations 
(4.21) s (4.24) end (4.25) are substituted in equation (4.16), 
Table 4.1 lists the final expression of torque obtained after 
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simplification for various corabinations of 'n' and *m' values. 
In the Table 4.1^ 'n’ corresponds to stator current harmonic 
and ’m’ the rotor cu3:ront harmonic when referred to the stator 
and p is a positive integer, i.e. 0,1, 2, 3 ... 

It can he seen from Table 4.1 that the torque spectrum 
has frequencies vvhich are raultiole of ’sBo)’. 

The frequency components of the type (3p+2) produce the 
flur rotating in opposite direction to the flux produced by the 
fundamental. The different torque expressions written in 
Table 4.1 can be v;ritten as one general expression if the 
frequency component of type (3p+2) and its phase is assigned 
a negative value. This implies that m or n takes the values as 


“2,3,4,- 




If this is done, the torque is given as 

T = 21 I M sin[(n--m)wt“(a''-cx" ) ] if (n-“m) is multiple of -3 

Q S I? C S x 

= 0 otherwise (4, on) 

where n and m takes the values as given by equation (4.27) and 


*>c 

c^s “ '-*^3 if n is positive 

= if n is iiegative (4.2,*t.) 

c:^ ~ ^-P if ra is positive 


if m is negative 


(4.1P,) 
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Thus if n and m are represented by values given by equation 

( 4 . 27 ) 5 is seen that the torque haraonic given by the 

difference of the stator and rotor current frequencies is 
produced. The equation for the torque (4,28) can be seen to 
be independent of the choice of the reference frame. 

4.2.2 Harmonic components of the torque for the induction 
motor fed by a cur.rent source inverter 

The stator current of the induction motor fed by a current 
source inverter is antisymmetric about mt = it:, i.e. i( 0 ) = 

-i(0 + -re). This implies the absence of the even harmonics 
in the stator current. The triplens are also absent because 
the sum of three phase current will be zero in this case. 

Thus the permissible values of stator harmonics from equation 

(4.27) are 

I 5 -5,7,-11,13,-17,19 (4.30) 

In a induction motor rotating at constant speed the induced 
rotor currents, referred to the stator will also have harraonics 
given by equation (4.30). Thus stator and rotor currents have 
odd and non triplens values only. 

It has been shown in Sec. 4.2,1 that torquo has contribu- 
tion from the harmonics vTnen (n-m) is triplen. Values of ’n* 
and ’m’ in this case, equation (4,30), are odd, and thus 
(n-m) is even. Thus the torque spectrum of the induction motor 
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■jeiao fed dy the current source inverter' have harmonics vdiich 
arc multiplies of six. 

In previous chapters j, the harmonics of the dc[ frame stator 
and the rotor currcents havo been obtained. The nth harr.onic 
an:’ nth hormonic of stator and rotor cl axis curraent can be 


i^l = Ij cos(n(ut - (4.31) 

^d? “ cos(ntot - ci^) ' (4.32) 


and on caao be related to i » c: and c:_ as follows. 
Substiturinp for the values for the values of current in equa-- 
tion (4. "7) from equations (4,1) to (4.3). 

i I^[cos(nnt“r:,) cosQ + 

C J. o S lo 

[cos[n((ot - §'”)-b 3 ) ] cos(9 - f”) ] 

[cos[n(wt - ) ] cos(e ~ ^-0] 

^ O s O 

Simplifying thiSj we obtain 

i^’i = Ig cos(no)t “ ci^ + (4.33) 

V'hero upper sign is v;hen n is of type (6p“l) and lo\';er 
o ';he r’ ifi s e . 

Co'.iparing equations (4.33) and (4.31) gives 



4.12 




a = a, + 0 

s 1 “ 


Similarly it can she shown that 


T = T 
-'-■r "-2 


i ^ i 


= cCq + 9 


(4.34) 
{A .35) 

( 4.36) 

(4.37) 
(4.37b) 




bstitating fron equations (4.34) to (4,37) in equation (4.28) 


sin[ (n“ra)a)t - (aj_~a 2 )] (4.38) 

Here n and ra takes value as given in equation (4.31). 

The signs of oc^ and (X 2 are reversed if n or m respectively has 
a negative value. 


It should be further noted that the harmonics of the 
pseudo rotor currents, and X 29 rather than rotor currents 
have been computed in the previous chapters. The dependence 
of torque on rotor currents can be changed to pseudo rotor 
currents by aquations (2.42) and (2.43). Substituting these 
in equation (4.14) gives 


T„ = Ih 




22 


■dl^ 


-dl 


X, 


M 

L 


22 


■q1 


)] 
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This impliGs that 


T 




X,'> 



(4.39) 


It is seen that this equation (4.39) is similar to the equa- 
tion (•'’;.!) except for tho change of i^2 ^q2 ^2 ^1* 

res oec tively . and X-j can be seen as ’d’ and q axis compo- 

nent froiTi equations (2.42) and (2.43). Thus a equation similar 
to equation (4.30) can be obtained in terms of the harmonic 
componeiics of the eseudo rotor currents as 


'T' 

i 


q 


'in' is 
values 


= 2I.,X^T;^ sin [(n“;.)o)t“(cij^'-c:^)] 

is mth harnonic ma>-_,nitude of X^, is 

positive aiid. negative of phase otherwise 
as given in equation (4.31). 

THODS FOR TFE COilPUTATIOI4 OF TIE TOR'qUE 


(4.40) 


its phase if 
’m’ take 


The expression relating the dq stator and rotor currents 
to the instantaneous electromagnetic torque is given by equa- 
tion (‘<5.1). In previous chapters, tv/o methods for the computa- 
tion of the rotor currents from a given stator current have been 
doscribod. 

The first racthod results in the time domain solution of the 
rotor current and the other gives the harmonic components of the 
rotor cue rent. 
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A method to compute the toro;ue produced by the squirrel 
cage induction is to substitute this time do.aain solution of 
the rotor current obtained, in equation (4.1) to get the time 
domain torque expression. This method is described in 
Sec. A .3.1. 

Other method uses the values of the harmonic components of 
stator and rotor current obtained. Equation (4,39) gives the 
expression for the harmonic torque produced due to the interac- 
tion of the nth harmonic of the stator and mth harmonic of the 
rotor current. For the system of motor being fed by the 
current source inverter, m or n takes values as given by the 
equation (4.31). For comoutation of particular ^torque harmonic 
Component, which is multiple of six (Sec. 4.2.2), all the 
possible combinations of the permissible values of *m’ and 'n‘ 
are taken which contribute to this torque harmonic. Their 
contributions to the torque can be obtained from equation 
(4.38), In Section 4.3.2, the procedure to evaluate all the 
possible combinations for a particular harmonic has been men- 
tioned 4 

4.3.1 Time domain solution of the torque 

In previous chapters, the method for the time domain 
solution of rotor currents has been given. There the cases of 
stationary and rotating rotors have boen dealt with separately. 
For the case of stationary rotor, the rotor currents expre- 
ssions have been obtained. This can be substituted in 
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equation (4.14) to obtain the torque produced for the stationary 
rotor case. 

For the case of rotating rotor the solutions obtained for 
and X 2 can be substituted in equation (4.40) to obtain the 
instantaneous torque values. 

It has been shown in Section 4.2.2 that the torque has 
harmonics to be multiple of six times the inverter frequency. 

It implies that the torque v/aveform repeats after every 60° 
interval. Thus the computation for this section has to be 

done only for any one of six intervals of the inverter. 

Frei^uency ' ‘solution "tKe 

The value of a particular torque harmonic can be computed 

through equation (4.40), For this computation, first of all, 

the possible corabination of tho current harmonics of the stator 

and rotor have to be obtained which contribute to this torque 

harmonic . 

The stator and rotor currents have infinite number of the 
harmonic components. The magnitudes of these harmonic go on 
decreasing as the number of harmonic increases. Thus though 
there are infinite number of the combination which will give 
the particular torque harmonic, the contribution because of 

I 

higher order current harmonics will be very insignificant. 

In the section the evaluations of stator and pseudo rotor 
currents harmonic combinations is done in a special sequence. 
This sequence has a property that at each step all the 
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permissible coubination of harmonic currents up to a specific 
frequency are coraputed. The torque value is updated with 
every new combination taken till it is observed that the 
contribution because of coming new pairs to the torque is 
less than a limit (say 0,001'/,) of the last updated value. 

This calculation of the torque harmonic has the following 
four steps. 

(i) To got the initial pair of stator and rotor frequency to 
start the torque computation. This has been shown in Fig, 4,1. 

(ii) To compute the torque value by a given pair of frequencies. 
This computation has been s'aovm in Fig. 4.2. 

(iii) To check when to terminate the present process of compu- 
tation, This is done by checking the contribution of the last 
pair to the torque; Further calculations end if this contri- 
bution is less than a limit (say 0,001><) and the two possible 
combinations with a component of pair 1 had been considered. 

If this is not the case we move to next step. 

(iv) This step computes the next permissible pair which con- 
tributes to the chosen torque harmonic. 

The flow chart of the procedure for computation of 1±ie 
chosen torque harmonic has been given in Fig. 4.3, 
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4.4 CALCULATIONS FOR A KNO/m MACHINE PARAMETERS 

A comput^x program is <i^v«lqp«d to cos^^ute torque harmonic 
using the flow chart of Tig* 4,3. The time domain torque is 
obtained via Sec. 4.3,1 approach. Tho listing of' -co^nputer pro- 
gram is given in the Appendix B. 

The average and sixth harmonic torque of induction motor 
produced by the stator current of Fig, 2.1 has been plotted in 
Fig, 4.4 as the function of rotor speed. It is seen that 
sixth harmonic torque is maximum at synchronous sp>eed. 

Fig. 4.5 gives the relative contents of various torque 
hamonics at throe rotor speed. It can be noted from this 
f igure that as the rotor speed is increased the torque harmonic 
components become mar e and more dominant . 

Fig, 4.6 plots the time domain torque waveform. This wave- 
form has been dra'wn for 60*^ interval only because the torque 
waveform repeats after every 60° interval (Sec* 4*2). 

Table 4,2 compares the values of torque harmonics obtained 
from following methods s 

(a) Using Sec, 4,3 method of obtaining torque harmonics from 
harmonics of currents. 

(b) Obtaining the harmonic components of the time domain 
solution of the torque* obtained through Sec, 4,2 approach. 
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(c) Obtaining tho valuo of torque harmonic by considering 
tho dominant current harmonics only. For tho torquo laarmonic 
froquoncy of 6 nw, the dominant current harmonics have fro— 
quoncy {6n-l)w and (6n+l)co because those interact with the 
fundamental to produce the harmonic torque. 

It is seen from Table 4.1 that all the three methods give 


tho similar results. 
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START 


Divide the nth torque harmonic component 
value into two halves, n/2 and n/2 




i Go on incrementing one value and decre- i 
I menting the other till neither is even or' 
j triplen i 
1 i 


The number corresponding to type (6p-l) 
is assigned the negative value 



The two numbers form the initial pair of ' 
frequencies ' j 

I 


V 


END :• 


Fig. 4 .1 s Step_l Computation of t he in itTa’* 
s.H?ti^_^air,ofXrj>^uen_ci7Bf^ 
and rotor current 
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start' 


Assign a boolean variable (say DONEKEV)] 
a value true ! 




j SuopcsG the input pair is (n,m) computes 
; nth harraonic of stator and mth harmonic of 
; pseudo rotor current 


j Compute torque value using equation (4.41) j 


I Add this computed torque value to the previous | 
j value keeping care of phase as this is vector i 
i addition i 

i i 



Toggle variable DONEREV 


i 

I 


Reverse the 
pair to 
(m,n) 


Pis » 4.»2 » Step 2 U pdating of torque with the 

nev / p a ir v alue 


No 




X 


- ‘MS DONEREV> 


True^ 

I 




( 


Yes 


■ END , 
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START ; 


(i) 


A boolean variable DONEPLUS 
is assigned true if chosen 
torque freejuency divided by 
six is odd else is assigned 
pulse 


(ii) DIFVAR 


- 6 


Step (i) of 
; Fig . 4.1 


I Increment 1 
1 DIFVAR by 1 
i number 6 | 

\ ^ ^ 

i j 


Increment j ! 


Decxement 

pair numbers i ■ 



pair numbers 

by DIFVAR i 




br DIFVAR 

' J ' 


Step (ii)' 
j of ’ Fig. 4. 2 


A., 


No 



IS \ 
DOTE PLUS " 
true? 


Yes 


\ 


\ 


/step (iii) >- 

/ 

/ 


Toggle 
DONE PLUS 


\ 


. continue 


Terminate 


'■ END 


Fig, 4 . 3 Fl ow chart for the computation of a 
p a rticular to rqu e harm onic 
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CHAPTER 5 

REDUCTION OF TORQUE HARivlONICS BY MODULATION OF D.C. 

INPUT CURRENT 


5 . 1 INTRODUCTION 

The induction motor being fed by a current source 
inverter produces harmonic torques. The harmonic torques 
can be reduced by modification of the stator current waveform 
This can be done either by modulating the dc current input of 
the inverter or by modulation v^ithin the inverter. In this 
chapter the former technique of modulation has been studied. 

The latter technique is referred as pulse width modulation in 
the literature [lO], 

The aim of the modulation is fco reduce the harmonic 
torques for the same average torque. The interest is directed 
more tov;ard reducing the dominant torque harmonics, 6th and 
12th. The higher order harmonics apart from being lower in 
magnitude are not reflected in the shaft speed because of the : 
large inertia of the mechanical system. 

The reduction in the torque harmonic can be noted by 
the increase in factor after modulation, where 

and 'T ' refer to average and nth harmonic torques respecti™ ' 
vely. This is so because torque can be seen to be proportionali 
to the square of stator current, from equation (4,38) . Thus 
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the inverter input current v/hen multiplied by the factor 
t(T^)no mod/( T^)mod] gives the average torque after 
modulation equal to case without modulation. This allows 
us to compare the relative contents of harmonic torques 
01 two torque spectrum of different average torques. 

It has been sho'wn that the modulating waveform has to 
fulfill a general requirement so as to produce a balanced 
three phase stator current. This modulating current should 
repeat after every 60^ interval, that is, the modulating 
current Jrrequency should be six times the inverter frequency. 

In this chapter the cases of modulating current vi^aveform 
of exponential or of consinusoidal nature have been studied. 
Having assumed a modulating current, satisfying the above 
requirements, the torque produced by the motor has been com- 
puted both in time domain and frequency domain (Sec. 4.3). 

Time domain solution is obtained by first evaluating 
pseudo rotor currents using method in Section 2.2. These 
expression when substituted in equation (4.39) gives the 
tim.e domain torque expression. 

For computation of the torque spectrum in frequency 
domain, the harmonic components of the stator dq current is 
obtained through a computer program. The harmonic components 
of pseudo rotor currents arc evaluated by proceeding as in 
Section 2.3. Then procedure in Sec, 4.3.2 can be used to 
compute the torque spectrum. 
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In Soction 5.2 the origin of the general requirement on 
the modulating waveform has been shown. In Section 5.3 the 
case of exponential waveform as the modulating waveform has 
been studied. In this case the amplitude and the time constant 
of the exponential waveforms are the parameters for the study. 

Section 5.4 the case of consinusoidal modulating waveform 
has been studied. In this case, the magnitude, frequency and 
the phase of the modulating waveform act as parameters for 
study. 

It should be noted here that in this chapter, for the 
purpose of simplification of analysis, the current source 
inverter has been assumed to be ideal. 

5.2 GENERAL REQUIREMENTS OF MODULATING V-/AVEFORJvl 

The current source inverter produces balanced three phase 
line currents. These currents have following property : 

(i) They are displaced from one another by 120° 

(ii) Each of the line current is antisymmetric (i.e., 

i(0) = -i(7t + e)) . 

Let the current during interval I in phase ’a’, be a general 
function of time f(t), where t refers to the time with origin 
at the start of an interval under consideration. For example, 
f(t) may be a exponential waveform as shown in Fig. 5.1, Then 
the intervals III and V of phases ’b’ and 'c’ respectively will 
also have the same current f(t ) due to property (i) above. 
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Interval IV phase 'a' current is - f(t) due to property 
(ii) above. Thus phase 'b‘ and * c’ have currents of type 
“ f(t) during intervals VI and II respectively. 

Thus during all six intervals , the conducting phases 
carry the sarae function of current f(t). Therefore} inver- 
ter input current should be same during all intervals to 
satisfy the above condition. This implies that the modulating 
current should repeat after every 60° interval, that is, 
it should have a frequency equal to six times the inverter 
frequency. 

5.3 EXPOIffiMTIAL MODULATION 

The general expression for the exponentially modulated 
inverter input dvc current, i(t), during a 60° interval of 
inverter can be written as 

i(t) = I + K(e^'^ - 1) (5.1) 

whore ’K’ and ’P' are the magnitude and inverse of time con- 
stant of the modulating waveform. 'I* is the unmodulated input 
d.c. current. Here the origin of time, t, is at the start of 
each of the intervals. The three phase currents obtained by 
this modulated inverter input d.c. current are drawn in , 

Fig, 5.1. To study the effect of the modulation on torque 
harmonics, both time and freciuency domain solutions for the 
torque are obtained. 



From the time domain solution of the torque it has 
been shown in this section that the torque harmonics can 
be reduced to zero for the case of stationary rotor, if 
K and jS are chosen to be ‘I* and (eqn» (2,18)) res- 
pectively. For the rotating rotor case it has been proven 
that it is not possible to get a constant average torque by 
using the exponential modulation. 

The effects of variations in K and p on the harmonic 
torque spectrum have been studied through frequency domain 
torque computation technique. The results obtain show 
that any exponential modulation improves the performance index, 
T^^/Tj^, for all rotor speeds below the synchronous speed. But 
in case of speeds greater then synchronous speed, there is a 
deterioration in this performance index if due to any exponen- 
tial modulation. 

5.3.1 Time domain analysis for the exponential modulation 

In this section, the analysis for the stationary and 
rotating rotor cases has been done separately. It has been 
shown that a constant average torque can be obtained for the 
case of stationary rotor but not dyring rotating rotor by 
exponential modulation, 

5.3. 1.1 Case of stationary rotor 

During interval I, Fig, 5.1, the currents are 
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where 

These 


= K e 
a 


gt ^ I* 


Ik = -(K + I*) 


"c =° 




I = (I - K) 


can be transformed to dq frame as 

idl = [K + I*] 

i,i = - r 

f3 


(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 


Tho time domain solution for the torque during this interval I 
can be obtained using equation (4.40). For this Xj^(t) and 
X 2 (t) have to bo calculated. 

For tho case of stationary rotor, i.e., co^ = 0, equation 
(2.51) gives 

P^2 = ^ ^iql 

This implies 


(a+p) X2 - ^ ai^ 

Substituting for i^-j. this equation from equation (5.6), 
gi. VQs 

(a+p) Xo = ei[K + I ] ( 5 . 8 ) 

^ ^22 
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Tho solution of this equation (5.8) is 

XaCt) = 

vvhero 

Kj K 

Ai = K?r 

■S€- 


(5.9) 


(5.10) 

(5.11) 


' (5.12) 

^ ^22 

and is a constant. 

Similarly, from equation (2.50) we obtain 

X^(t) = C 2 (5.13) 

whore 

1 ^ ^ 

Di = - p • ri+F) 

j, *^1* 

f3 ^ 

and C 2 is a constant. 

Substituting for i,^j_> ^ql* ^1 ^2 equations (5*6), 

(5,7), (5.14) and (5.9) into equation (4.40), wo obtain tho 
torque Gxprossion during interval I, after simplifying as 


(5.14) 

(5.15) 
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The time dependence of torque in above equation can be removed 
if K and p are chosen as 

K = I } p = a (5.16) 


For this choice. 


T 


q 





I 


(5.17) 


that is, a constant torque is obtained. 


The values of and C 2 depend upon the initial conditions 
and can be computed by the analysis, as shov«i in Sec. 2.2. 
Using the method of Section 2.2, the follov;ing values can be 
obtained 


'1 1 


C2 I 


f3/2 


(4^3) I 


i -ir3/2 
L 


(k2“4) 


L f3 




x]; = (| - 


(5.18) 







+ B 


1 


^3 " -^1 


(5.19) 
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Table 5.1 lists the expressions of ^ql^ ^1* ^2 

T with values of K and p from equation (5.16) for intervals 
II and III, Here C^, and are constants which depend 
on the initial conditions and can be calculated as in 

Section 2.2, In this, 62 and t^ are as defined in Section 2.2, 

The torque is seen to be constant for intervals II and III for 

K,p choice of equation (5,16). Intervals IV, V and VI are 

similar to intervals I, II and III respectively because of 
symmetry of current about wt = -n:. 

Substituting, from equation (5.16), in equation (5,1), 

gives 

i(t) = I (5.20) 

This equation (5,20) gives the form of the input current 
during an interval, which will produce a constant torque for 
the stationary rotor case. 

5. 3. 1.2 Case of rotating rotor 

The pseudo rotor currents are related to stator dq 
currents from equations (2.52) and (2.53) as 

(p2+2ap+a2«2)x^ = (5,21) 

(p +2ap+a ~ ^l^^dl ”” ^l^r^ql ^l^^dl (5.22) 

Let us confute the torque expression for a period during inter- 
val I i.e. for *u)t' values between 0 and 7t/3. 
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Using values of this period in (5,21) 
and (5.22) from equations (5,6) and (5,7), the expression 
for and X 2 are X 2 are obtained as 


^2 
where 


= e cos(ti>^t + 0^) 

Sis 

+ ^ B 2 

(5.23) 

= Cg e”^"^ cos(a)^t + 0^) 

+ Ag e^'*' + Bg 

(5.24) 


K Qi 

2 +2ap+a^4-o)^ 


(5.25) 


= ^l^“r *" 


B 


qo 




igjtsaeagj!? ^ 

^ (3^-K»p 


qj = - p) 



( § +2a^+a +b)p 


0) 


^3 ~ ^l^fS + p) 


B 


q4 




3 

^4 


K^(f§ + a) 


(5.26) 

(5.27) 

(5.28) 

(5.29) 

(5.30) 

(5.31) 

(5.32) 
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equations (5.23) and (5.24) are constns 
which depend on the initial conditions and can be computed 
by the analysis as done in Section 2.2. 

The expressions for and are substituted 

fro» equations (5.6), (5.7), (5.23) and (5.24) into equation 

(4.40) to get the expression for torque. After slrapllf ication 
this equation is given as 

Tq = e^^ Cg cos(m^t+ 03 ) ^4 + 

+ (KB4+A4I )e^VCgI e cos(a)^t-f-03)+B4l*] (5.33) 


v;here 




64 = 3 " 


4 ^l^r^ 




(a-^) 


and Cgf03 are obtained from relation 


(5.34) 


(5.35) 


c 

Cg cos(a)^t+03) = Cj cos(w^t+03) cos(u)^t+02) (5.36) 

That iSy Cg and 0^ depend upon the initial conditions. : 

To remove the time dependence of equation (5.33) we have the 

choice of two parameters K and p. Sinco there is a term of 
23 1 

type , to make equation (5.33) time independent, it’s 
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coefficient should be made zero by choice of K or p. This 
is possible only either K or is zero. Since we are 
Dealing with rotating rotor case, ^ 0, For exponential 
modulation K ^, 0 . Thus (K.A4), the coefficient of e^^^ cannot 
be reduced to zero. 

Thus it is not possible to get an exactly constant torque 
for the rotating rotor case, using the exponentially modulated 
waveform. 

Though it is not possible to make equation (5.33) inde- 
pendent of time v/ith K and p as parameters, but it is evident 
from this equation that certain choices of K and ,8 may give 
lower time variations in torque, then others. This can be 
seen in a better v/ay through torque spectrum in frequency 
domain. This has been studied in detail in next section 
through the frequency domain analysis of the torque spectrum 
duo to the current waveform of Figure 5*1. 

5.3.2 Frequency domain analysis for the exponential modulation 

The effects on chc torque spectrum duo to the variations 
in the exponentially modulated waveform, Fig. 5.1, have been 
studied in this section. The torque spectrum produced by this 
stator current waveform has been computed using the method of 
3ec. A. 3. 2. 

It is seen from equation (5,1), that for exponential 
modulating waveform there are two control parameters, p and K. 
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In this section the effects on torque spectrum due to the 
variation in the values of p and K of the modulating waveform 
has been studied for various values of rotor speeds. When the 
effects due to variations in ’p’ are studied, ’K’ is assigned 
a value equal to ’I*, which gives no torque harmonics for 
stationary rotor case, (Sec. 5.3.1.1). Similarly, when effects 
on torque because of variations in ’K’ is studied, ’p* is 
assigned the value 'a’, as this value of 'P* removes torque 
harmonics for stationary rotor case. Sec. 5. 3. 1*1. 

5. 3.2.1 Analysis with the variations in the time constant of 
the modulating exponential waveform 

To study the effect on the torque spectrum due to variation 
of p, the inverse of modulating waveform time constant, the 
amplitude of this waveform, K, is taken as I. It is noted that 
this value of K and p equal to ’a' give a constant torque for 
stationary rotor case. For this choice of K, the expression for 
the modulated inverter input current is given by, from equation 
(5.1) as 

i(t)=Ie^^ (5.21) 

For the current the torque spectrumccan be obtained for 
various values of p , at various rotor speeds. Suppose we are 
interested in reduction of the 6th harmonic of the torque. 

Fig. 5.2 plots the performance index (Tay/l^) function of the 
rotor speed. refers to 6th harmonic torque and to the 
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average torque. In Fig. 5,2, p, the inverse of the time 
constant of the modulating waveform has been taken as para- 
metier. The case p = o, in Fig, 5.2 corresponds to the no 
modulation case. It is evident from this figure that 

(i) For the rotor speeds less then the synchronous speed 

of the motor there is improvement in the performance with aaf* 
value of of the modulating waveform as compared to the no 
modulation case. As against this, for the rotor speeds 
greater than the synchronous speed of the motor there is no 
improvement with the exponential modulation, rather there is 
deterioration in the performance. 

(ii) It is seen from Fig. 5.2 that there is a optimum value of 
p , at particular value of rotor speed which gives the best 
performance. For rotor speeds less then synchronous speeds, 
this optimum value of p increases as the rotor speed increases. 

(iii) (a) The performance index is a sharp function 

of p at rotor- speeds much loss, then synchronous speed. Howeverj 
at higher speeds, that is, closer to synchronous speed (T^^/T^) 
is not very sensitive to values of p. For example j at rotor 
speed of 300 rpm, (T^^y/T^) for p = a is 20.3 and for (p=2a) is 
31.4, Siinilarly at 200 rpm when p changes from a to O , 

reduces from 36,0 to 17,4, At higher speed say at 
560 rpm, (^.^y/^^) for P = 12a and 16a is 5,87 and 5.06 res- 
pectively. At 580 rpm, (T^^/T^) changes from 1,41 to 1.98 
for change of p from 12a to 16a. 
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( b) In view of above, it is possible to keep p constant in 
the normal operation region of induction motor, i.e. between 
approximately 90/. to lOO;:^ of the synchronous speed and still 
obtain near optimum value of 

(iv) As the rotor speed tends to zero the optimum value of 
’p’ approaches the value *a', as obtained in Sec. 5.3.1, 

5. 3, 2. 2 Analysis with the variations in the amplitude of 
the modulating exponential -waveform 

To study the effects on torque spectrum due to tho varia**- 
tions in the amplitude of exponential modulation, K, the value 
of p is kept constant. This constant value is taken as ’a*, 
as this choice of p gives constant torque for stationary rotor 
(Soc. 5. 3, 1,1). In this section, the interest has been focussec 
on the reduction of the sixth harmonic torque with variation 
in K. 

Fig. 5.3 plots the performance index (Tgy/Tg) as the 
function of rotor speed, varying from zero to synchronous 
speeds. The case K = 0 in the figure corresponds to the case 
with no modulation. It is evident from this figure that the 
effects on torque spectrum due to variations in K are similar 
to ■ the variations in torque spectrum v^ith p as a parameter. 

That is, for every value of rotor speed, there is an optimum 
choice of the value of K, which gives us the maximum perfor- 
mance index, (T .JT^) . Similarly, this index is a strong 

O V O 
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function of ’K* at rotor speeds much less then synchronous 
speeds. Hov/ovor at higher speeds, that is, speeds closes to 
the synchronous speeds is not a very sensitive to 

the values of K. In view of this, it is possible to keep K 
constant in the normal operation region of the induction motor 
and still obtain a near optimum value of 

It can be seen' from Fig. 5.2 that as the rotor speed 
increases, the optimum value of K for the best performance 
increases. For rotor speed approaching zero, this optimum 
value of K approaches the value I, invertor input dc current wit 
no modulation. This is in fact the optimum value of K, which 
is obtained for the case of stationary rotor (Sec. 5. 3, 1.1). 

5. 3. 2. 3 Optimum pair of (K,p) for a given speed 

It has been observed in Sec. 5. 3. 2,1 that for a given rotor 
speed and a assumed constant value of K, there is a value of p, 
•which gives the maximum performance index Similarly, 

Sec. 5. 3. 2. 2 shows that v^hen p is taken as constant, this opti- 
mum property is observed with respect to the value of K as well. 
It has boon shown in this section that there is a optimum 
choice of the pair (K,p) v\rhich gives a maximum performance index 
at a speed. 

To observe this , die maximum T^^/T^ value is obtained for 
various choices of p by varying K. The plots of 
as a function of p and (K) optimum as a function of p has been 
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plotted in Fig, 5,2 for tho case of rotor speed = 500 ipn. 

F-rom those plots it can be observed that, 

(i) There is an optimum value of the pair {K,p) which gives 

tho largost value. For this pair, the 

is very largo and it can be assumed that a near constant tor- 
que is obtained with this choice of (K,?) pair. 

(ii) As the value of p increases, tho value of K which gives 
tho maximum value of (T^y/X^) for tho chosen p , goes on 
decreasing. 

Fig. 5.3 plots the optimam (K,p) values as the 
function of speed. From this plot it is soon that as tho rotor 
speed increase the optimum value of p goes on increasing. It • 
is equal to 'a* , at zero rotor speed, 

5.4 COSINUSOIDAL MODULATION 

Tho general expression for the cosinusoidally modu- 
lated inverter input dc current, i(t), during a 60° interval 
of invertor can bo written as 

i(t) = I + Am cos(to t + a ) (5,38) 

whore A^, (x^ are the amplitude, frequency and tho phase 

angle of the modulating waveform. ’I’ refers to unmodulated 
d.c. current. The effects on tho torque harmonics duo to such 
a modulation has been studied in this section with both the 
time and frequency domain torque solutions. 



5,18 


From the time domain solu'fcion of the torque it has shown 
that it is not possible to get a constant average torque with 
a cosinusoidal modulation. The effects of variations in 
and on the harmonic torque spectrum have been studied 
through the frequency domain torque computation technique. 

5.4.1 Time domain analysis for the cosinusoidal modulation 

In this section, it has been argued that it is not possibi< 
to obtain a constant torque v;ith cosinusoidal modulation. The 
cases of stationary and rotating rotors are dealt together. 

Lot us compute the torque expression for a period of 
interval I i.o, for ’wt’ values given by 

o < wt < ^ 

For this period, the three phase currents for the dc input 
current given by equation (5,38), are 



(5.39) 


(5.40) 

i =0 
c 

(5.41) 

These phase currents can be transformed 

equations (l,9)and(l*10) as 

to dq currents from 


(5.42) 


(5.43) 
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Tho expressions for and can be obtained by solving 

tions ( 5 , 21 ) and (5.22) v/ith above values of i , and i « 

ql q2 


gives 


oqua- 

This 


X3^(t) = cos(w^t+ 04 )+D^ cos(o)j.^t+a^)+D2 

(5.44) 

X2(t) = o"®'*' cos(w^tt0^)+D2 cos(Wj.,jt+a^)+D^ sin(w^t+a^)+E2 

(5.45) 

whore 



q2 


q 

D 


2 


D 


1 

3 


q4 


o o o 

= Am [q2(a‘"+m^+w^)-2aqu}j^]/E3 
= -[4a-ta^^+(a^+m^-w^)^] 

_ ^1 
- T3" 

= Am[2aq2Wjjj+q(a^-to^-w^)]/E3 

qg I 

= '' ~w 

a 

= Am[ q^i a^+w^-w^) -2aq^w^]/E3 
= K^{a+ ^f) 


(5.46) 

(5.47) 

(5.48) 

(5.49) 

(5.50) 

(5.51) 

(5.52) 

(5.53) 


95 “m 


(5.54) 
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= Aiii[2aq4Wj^+q5(a^-H>i^-w^)]/E3 
r- _ *^4 ^ 

a^+wj 


(5.55) 

(5.56) 


C-jQ and Cj^j. constants v\^iich depend on the initial con- 

ditions. These can bo evaluated by the analysis as done in 
Sec. 2,2. 

The expressions for and obtained are substi- 

tuted in equation (4.40) to get the expression for torque, Afte; 
simplification this equation is given as 


^q “ '^''^c^^l2 ^ cos(a)j,t+0^)+lD^ cos (a)^t-i*a)+IDg^ sin(u)jj^t+a) -f 

+ cos(oj^t+0^) cos(Wjjjt+a) -i'(^)M D5 cos(2a)jj^t+2a) + 

H-(^)jVl sin(2w^t + 2a) + (^)M (D^+D^)] (5.57) 


Vifhere 

(5.58) 

(5.59) 


C^o and 0^ are obtained from equation 

cos(a)^t+0^) = C^Q cos ((0^1+0^)+ cos(a)^t+0^) 


D5 = 4M Kj_ Wj,(a -»^r“Wj^)/(3ETL) 
Dfi = 8S a Kj_ u>^ V(3 Ei) 


(5.60) 
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If equation (5.57) is to give a constant value, we should 
satisfy follov./ing relations 

I[D5 cos(u)^t+a) + sin(a)j^t+a) ] = O (5.61) 

M[D^ cos(2{j>^t+2a) -i- sin(2w^^t+2a) ] =0 (5,62) 

!/iCio = 0 (5.63) 

ICj^2 ^ (5.64) 

equation (5.62) can be- vsfritten as 

M''/'(D^+D^) cos(2Wjj^t+2a+0^) ==0 (5.65) 

where 

tan 0j = -(D^/D^) (5.67) 

Equation (5,65) is satisfied only if M = 0 or and are 
zero (5,68) 

M = 0 implies the case of no modulation and is not of interest. 
Thus , wo should have 

— 0 (5.70) 

=0 (5.71) 

and from equation (5,63) 

^12 


= 0 


(5.72) 
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Equations (o«70) ato (5,72) arc the only possible condition 
which satisfy equation (5.61) to (5«64), for cosinusoidal 
modulation. Substituting those in equation (5,57) gives 

Tq = 0 (5,73) 

This corresponds to the trivial, case of no stator current. 

Thus it is not possible to roduco current harmonics to zero, 
using the cosinusoidal modulation. 

5.4.2 Frequency doraain analysis for the cosinusoidal 
modulation 

The effects on the- torque spectrum due to the variations 
in the amplitude, frequency and the phase of the cosinusoidal 
modulating waveform have boon studied in this section. The 
torque spectrum produced by the assumed stator current is 
computed using the method of Section 4.3.2* 

The torque spectrum for the cases with of equation 

(5,38), equal to 3,6,12 and 18 times tho inverter frequency have 
been studied. When s tudying tho nature of torque spectrum it 
is found that for a particular choice and there is a value 
of which gives tho maximum performance index (T^^/T^) . Also 
thoro is a optimum combination of amplitude and a frequency 
which gives the lov^est torque harmonics. 


The results have boon summarised as follov/s s 
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~ ^ ° TablG 5,2 gives the 6th and the avorago 

torque for various araplitudos and phases for the case of a)„ 
equal to 6to and a fixed constant rotor speed. This shows that 

(a) There is a distinct phase angle for a every choice of 
amplitude which gives the maximum performance index for 6th 
harmonic . 

(b) There is a optimum amplitude which for a particular phase 
angle gives zero sixth harmonic torque. 

(c) As the amplitude of modulating waveform increases the 
phase angle at which the minimum 6th harmonic torque is obtained 
decreases . 

(d) The variations in ‘the average torque produced duo to 
modulation are not much significant. 

(ii) In Table 5,3, the effects on 12th harmonic torque for the 
case of 0 ) = 6a} has boon taken up. It can be noted that wq arc 
able to obtain the reduction in the 12th harmonic torque with 
modulating current frequency equal to six times the invertor 
frequency. In this case also the property of the optimum phase 
and amplitude has been observed. 

(iii) = 30} 3 Table 5.4 lists the results for equal to 
three times the invertor frequency. This typo of harmonic can 
bo present in the inverter input dc current due to imperfect 
filtering. It is noted that there is an improvement in the 
ratio of average to 6th harmonic torque and the feature of 
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optimum amplitudo and phase angle is present and is similax 
case (i). But for this case the change in the average torq\ 
is found to be significant. Even for a particular amplitude 
of the modulating waveform, there is a significant variation 
in the average torque with the variations in the phase angle 

(iv) = 12o) 5 The reduction in 6th harmonic torque by- 
modulating waveform with equal to twelve times the invert* 
frequency has also been studied. Table 5.5 lists the result 
obtained for this analysis. It is observed that in this case 
the variation in the sixth harmonic torque due to the varia- 
tions in tho amplitudo .and the phase of the modulating wavefo 
is small, as compared to the frequency of 6 uj. Also, hero a 
largo scale modulation is required for better performance. He 
no optimum property with respect to tho amplitude of modulati! 
current is observed in tho possible range of amplitudo of modi 
lation, 

(v) Table 5.6, lists the result for tho reduction of 6th har- 

monic torque with = 18o). Here again the variations in 6th 
harmonic torques are observed to bo small. These arc oven 
smaller then for the case (iv) above. Here also no optimum 
property is observed. i 

Thus it can bo said that in case the cosinusoidal modula- 
tion is being used, tho modulating waveform of the frequency 
six times the invertor frequency gives the best control for 
obtaining tho desired optimum performance. 
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Fig 5.3; Performance index Vs. rotor speed 
with fi- a and variable K 
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1.000 

72,000 

1 1,84115 

2,54987 

5.9 

3,000 

lOH.OOU 

15.26033 

1.03751 

9.3 

3.000 

1 41.000 

15.46141 

0,80943 

19.1 

3.000 

1 110,000 

15.50565 

0.14185 

109,3 

3,000 


15,46341 

0,63600 

24:i 

3,000 

252,000 

15,26667 

1,77166 

8,6 

3,000 

283,000 

14.85130 

3.01507 

4.9 

3.000 

32 1,000 

14.37406 

3.82979 

3.0 

4,000 

0,000 

14.07609 

4,54330 

3.1 

4,0 00 

36,000 

14.39424 

3:79217 

3.8 

4,000 - 

4,000 

72,000 

108,000 

15,08982 

15,649S9 

2,72842 

1:66756 

5,5 

9,4 

4.000 

4.000 

144.000 

180.000 


0:90505 

0:75962 

20.9 

4,000 

4,000 

216,000 

232,000 


0,86196 

i:9i39j 

10,4 

0,2 

4,000 

268,000 

15,10421 

3,4BOO? 

4:3 

4,000 

324,000 

14,40652 

4.51192 

3.2 



3 I'l^'f.r- 1 ;■/, Mf.3, jiT 

WITH iDCse.l) 


KTIH 
HUT{jp, 
?, 1 P« 


MOD. FFEO. 
SPFED=520, 


AS 

PPM 


Afi!> OF 
'■\fjV(AttP) 

<>,00{; 

0,‘3<)0 

0,500 

0,500 

0.500 

0.500 

0,500 

0.500 

O.SOo 

0,500 

0. 500 

1,000 
1 . 000 
t.ooo 
1 .000 

1,000 

1 .000 

1,000 

1,000 

1, ^000 

1.000 

2*550 

2.000 

2.000 

2.000 

2,000 

2,000 

2:000 

2.000 

2.000 

2,000 


-IF 

''lrjr)(nt:G) 


AVmhCE 
TOP 01 IF 


'.“'JC 14.64299 


0,0 00 

36.000 
72.tiO( 

I 08 ,( 0 >U 

1 44.000 

180 . o <‘0 

2lfc,oo<» 

252.000 

289.000 

324.000 

0,000 

36.000 
72,U(H) 

108.000 

144.000 
IwO.OOO 

216.000 

252.000 

288.000 

324.000 

. 0.000 

36.000 

72.000 

108.000 

144.000 
18O.UO0 

216.000 

252.000 

288.000 

324,000 


14.67053 
13.29620 
12,47106 
12.46395 
13,27705 
14. 64560 
16.07634 
16.99415 
17,00232 
16,0972? 

14*'^2843 

12.03599 

10,47836 
1C. 465 19 
11.99939 
14.67897 
17,59626 
19.52454 
19.54195 
17.63972 


14.93532 
9,77526 
“ 03065 
00855 
70892 
14,83639 
20.89581 
12301 
16206 


7, 

7, 

9, 


25 

25 


20.98959 


HAPFOiaC 

XfIFOUF 

4,08537 

5,66336 

4,93793 

3,946P1 

2.99474 

2.42174 

2.55542 

3.3372? 

4. <*1 237 
5.3741 7 
5,84672 

7,21695 

5. ?46?3 
3.92273 
2.15428 
0.H5713 
3 .48910 
3,01649 
4,97940 
6,83709 
7.70583 


10,18541 

7.01585 

3,93735 

1,58291 

1,99219 

3,71072 

4,68565 

6.87766 

10,24406 

11.71798 


Tav/Tnar 


3.6 

2.6 

2.7 

3.2 

4.2 

6 K 

5:7 

1*s 

3.9 

3.2 

2.8 

2,0 

2,1 

2:7 

4.9 
14.0 

9.9 

5.8 

3.9 

2.9 

2.3 

1,5 

1.4 
1.8 
4.1 

4.9 
4,0 

?•§ 

2*7 

2.5 
.1 .8 



^TABLE 5,5 


FOR REDuriDl !! 
1 ? Tinea InV , 
w X r N 


*F 6 HAH. 
EHEO. AT 
THC = 6 AMPS 


WITH 

RilTlJR 


A HP OF 
MUf)( A.MP) 

0,000 

0.500 

0,500 

0.500 

0.,500 

0.500 

0,500 

0.500 

0,500 

1*500 
t ,500 
.1,500 
1.500 
1.500 
1.500 

l-ioo 

1.500 


,500 
,500 
> 50(1 
300 
,500 
,500 
300 
300 


4.000 

4.000 

4.000 

4.000 

4.000 

4.000 

4.000 

4.000 

5.000 

5.000 

5.000 
5,000 
5,000 
5,000 

5.000 

6.000 
6,000 
6.000 
6,000 
6,000 
6,,#0 


Ai»GijE Ih 


average 

TORydE 


O.'JvO 14,6 42H9 


U, «>()() 

72 1 000 
100,000 
160.000 

252.000 

200.000 

324.000 

p.OOO 

36.000 

72.000 

1 00 , Coo 

180.000 

252.000 

288.000 

324.000 

0,000 

36.000 

100.000 
180,000 

252.000 

280.000 

324.000 

^P.OOO 

36,000 

'/ 2.000 

100.000 
100,000 

252.000 

288.000 

324.000 

36,000 

. 72,000 

100.000 

180,000 

252.000 

324.000 


1 

1 

14 

14 

14 

14 

14 

14 


, 61 425 
.62053 
.63702 
.65620 
.67650 
,65697 
.63776 
.62106 

37111 
30369 
.{>491 1 
.207^ 
,75707 
.70890 

.65152 

.59542 


1 1,54709 
1 4,59072 
14,69325 
H. 7901 1 

14.79303 

14,69754 

14.59406 

14,54690 

14.63111 

14,81956 

14.97523 

15.04492 

14,97920 

14,82707 

14.63753 

14,68793 
14,94383 
,l|897 

Iflll 


* ** 1 

If- 


6.000 


14,69685 

0.000 14.613S4 

7g.0pp 15,10014 
If, 13498 

isjol JS 06 S 57 

124,000 


Mun, FRKO. 
. SPEEn = 520 , 


HARMOnTC 

THROUE 


A 6 
PPM 


4..)0537 


1 * 12^23 

3.99079 
H , 09 b 40 
19314 
25021 
10306 
990 23 
92519 


3,o0045 
3 . 65725 
4,19937 
4,45090 
4.50383 
4,23337 
3,69725 
3,61043 

3,28704 

4;39372 

!- 2 I 433 

4,92205 

4.48632 

3.91029 

3.33770 

2,83720 

3.05196 

l•!2?53 

5,30255 

5,43660 

5,07838 

4.16614 

2.99266 

2,55568 

:?!??! 

I-IP’’ 


9462 





Tav/r lar 


3.6 

3.7 
3.7 
3,6 

3.5 

3.5 

3.6 

3.7 

3.7 

4.0 

3.0 
3.5 
3,3 
3.2 
3.5 

3.8 

4.0 


4 

3 

3 

3 

3 

3 

3 

4 

5 
3 
3 
2 
2 
2 

3 

4 


4 

fl 

3 

1 

0 

3 
8 

4 

1 

8 

1 

« 

8 

9 

6 

9 


5.7 

3.7 
2.9 

2.7 
2,6 
2,7 
5.2 

6.4 
3.6 
2*7 

2.5 

2.5 
5,4 


TAHLE 5.6 


FlJH 

18 


Pt'nuCTTi.K! 
times T‘,y 


i!F 6 H.\R. WITH MUD, KPEO 
. FKLQ. AT RUTUD SPEEn=52(> 
tDCse.O AMPS 


AS 

RPf’ 


Ahp or 

np 

average 

harmonic 

Tav/Tfiar 

MDIH AM?’) 

M'tncnbG) 

TUR-OUE 

TOROUE 


O.non 

O.OGO 

! 4.64299 

4,08537 

3.6 

0,50(5 

0,000 

14,62872 

4,02507 

3.6 

0,500 

JO.OuO 

14,63185 

4,04398 

3,6 

0,500 

72.000 

14.64013 

4,07950 

3,6 

0,500 

10 a. 000 

14,64991 

4,11604 

3.6 

0.500 

144,000 

14.65743 

4.14048 

3,5 

0,500 

1 60 .Ot'O 

1 4.66030 

4.14b0u 

3.5 

0.500 

216,000 

14.65771 

4.11130 

3.5 

0.500 

252.000 

14.65039 

4,09y9b 

3,6 

0,500 

288.000 

14.64(S65 

4,06181 

3,6 

0,500 

324.000 

14.63219 

4.03218 

■ 3.6 

1,500 

0,000 

14.60927 

3,90550 

3.7 

2,500 

0,00 0 

14,60194 

3,78751 

3.9 

3,500 

0,000 

14.6U675 

3.67132 

4.0 

5,000 

0,000 

14,63670 

3.50100 

4.2 

6,000 

0,000 

14.67183 

3,39059 

4.3 

6.000 

36.000 

14,74586 

3.91350 

3,8 

6,000 

72,000 

14,90824 

4.51016 

3.3 

6.000 

108,000 

15.02815 

4,81282 

3,1 

6.000 

144,000 

15,05698 

4.84931 

3.1 

6,000 

180,000 

15.05079 

4,82955 

3.1 

6,000 

216.000 

15.05619 

4189362 

3.1 

6.000 

252,000 

15.03139 

4.83267 

3,1 

6, you 
6.000 

288,000 

324,000 

14.91709 

14.75413 

4,33743 

3,50974 

3.4 

4.1 


6,1 


CliAPVdR 6 

CALCULATIOM OF CURRENT PROFILE FOR SPECIFIED 

TCRyUE V/A'/EFOKvl 


6.1 INlMODUCTIOi: 

In the previous chapters the form of the- s.tator current 
of the motor has been assumed and then the torque wavefomi 
has been obtained. In this chapters the inverse problem of 
the computation of the stator current to obtain a assumed 
torque waveform has been studied. This stator current is 
assumoo to bo produced by an ideal throe phase current source 
inverter whoso input dc current is being modulated, 

A nonlitooar second order equation has been derived for an 
interval of the invertor. The solution of this equation gives 
the stator curiont v/aveform for the assumed interval. Tho 
stator current waveform needs to bo computed only for one 
interval as the invertor input current needs to be modulated 
at six times tho inverter frequency (Sec. 5.2). Thus the 
s’oaror current during complete invertor cycle can bo obtained 
from tlr- seslution of stator current during any one interval, 

Tho equation to be solved is of tho second order, Tho 
solution of this equation requires two boundary conditions. 

The- so have been obtained, using the fact tliat under steady 

$ 

3 Ur to, xhe- rotor mmf must move by 60^ in space in each t/6 
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intorval so that each intorval is similar. 

In this chapter, the stator current \ 7 avoform roquirod to 
obtain a constant torque- has boon derived , For this case, 
the nonlinear second order equation reduces to two first 
order nonlinear equations. The analytical solution of 
those equations for the stationary rotor case has been 
obtained. This gives a exponential modulation as obtained 
in Section 5.3.1. For the rotating rotor case, those equa- 
tions have been solved through tv/o techniques, namely phase 
plane analysis and nuraorical analysis. The numerical solu- 
tion has been done by a computer program. The solutions for 
the stator current which give the constant torque have been 
thus obtain-.id for various values of the rotor speeds. 

In this chapter the interval of the invertor for which 
the analysis has been clone is intorval I (Fig. 2,1). It 
should again bo noted here that the inverter has boon assumed ■ 
to bo ideal for the analysis. 

6.2 DERIVATION OF EQUATIONS FOR STATOR CUR lENT TO PRODUCE 

DESIRED TORQUE '/AVEFORU. 

In this section, th^ equation for stator current wave- 
form during interval I is obtained, which gives an arbitary 
torque waveform T (t). This equation is a second order non- 
linear equations. The two boundary value conditions required 
to solve this equation have also been obtained in this section. 
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6,2.1 Derivations of nonlinear equation 

Lot us assunc- a General form i^(t) of tho ’a’ phase 

a 

current during interval I, v/hich gives the desired torque 

v;avofori;i T (t). During interval I. for the ideal inverter 

q 


“^a 


ic = 0 .... 

Thus from equations the dq stator currents can be 

obtained as 


( 6 . 1 ) 

( 6 . 2 ) 


^dl = 


"ql 


f3 


(6.3) 

(6.-^1) 

Substituting these values of i^^. i^^ eqv’ation (4.4^? 

vvG obtain 

X. 


V(t) = -H i (t) [X. + ™ 3 

q c a ^ Y3 


(6.5) 


In equation (6.5), T (t) is an assumed function of time and 

i (t) has to bo evaluated, 
a 

Fron equations (2.52) and (2.53) 


(o2+2ap+3^+<o2)x^ = K^al^j_ + 


(p2+2a?+a2+to2)X2 = ~ 


dl 


( 6 . 6 ) 

(6.7) 
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where is rotor speed and a, are given by equa 
ana (2.53a) respectively, 

r ?. oia equa Lions (6.6) and (6.7), v/e have 
(p2+2ap+a2+u2)(x^ + p) = fr' + 


tion (2.18) 


rr> " Vr [idi - 


Substituting for i^^ and i^-j in above equation fron equations 
(6.3) and (6.4), we have 


(p^-+2ap+aW)(X^+ p) = I Kl'^r^a 


( 6 . 8 ) 


Substituting for i_(t) in above equation from equation (6.5) 

3 . 

gives 


{p^+2ap+a^+w^)(Xj^+ p) 


*l“r 

ralip' ‘■'X, 


(6.9) 


I'his can be written as 


( p^4-2ap+a^-HA)^) (y(t) ) 


■’yTtT 


(6.10) 


where y( t) = (X^ + j^) 


( 6 . 11 ) 


Equation (6.5) using equation (6.11) gives 
i^(t) = ~T^(t)/(i\ y(t)) 


( 6 . 12 ) 
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Equation (6,10) can be solved to obtain y(t). This on 
substitution in (6.12) gives the form of stator current 
which gives the desired torque waveform T (t). Equation (6,10 
being of order two requires two boundary value conditions. 
These have been obtained in Sec. 6.2,2, 

6.2,2 Derivation of the boundary conditions 

The values of dq currents at the end of the inverter 
period are related to the initial values by the equation 
(2.77) as 


X 2 (T^) * (1/2) X^io) - (^3/2) X^(o) (6.12) 

Xi(Tc) = (f3/2)X2(o) + (1/2) X^(o) (6.13) 

where T^ is the time of interval I, that is, (u/3co), co being 
the inverter frequency. Thus, Xj_(T^) and X 2 (T^) correspond 
to the pseudo rotor currents at the end of interval I and 
Xj^(o) and X 2 (o) at the start of the interval. Those relations 
are used to compute two boundary value conditions in terms of 
y(t). 

From equations (2,50) and (2.51) 

X^(o) =6)^ X2(o) - aX3^(o) + iqi(o) 

^2(0) "^r ^ 1 ^°^ “ 0X2(0) + 


(6.14) 

(6.15) 
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Here 


XAo) 


(6.16] 


Xr)( O) — 


?roQ equations (6.14) and (6.15) v/e have 


(6.17) 


( ( o) + 


Xo(o) 


X,(o) 




“Y*;^”) ~ - “Y^””) *'a(Xj^( o) + 


^^•1(0) 

Using definition of y(t) from equation (6.11) and equations 
(6.3) and (6.4) in above, we obtain 


y(o) = a)^(X2(o) 


Xi(f>) 


■)-ay(o) 


(6.18) 


y(o) = 


(6,19) 


From equation (6.12) and (6.13) 

X..(T^) , Xo(o) X, (o) 

(Xi(Tc)+ -^-—0 = 2 CXjl(o) + -p— ) + p (^(o) - -p^) 


Using the definition of y(t) from equation (6.11) and substi- 

X,(o) 

tuting for (X2(o) - froBi equation (6.18) in above we 

obtain 


y(T ) = 1/2(1+ ^^) y(o) + ^ y(o) 


( 6 . 20 ) 
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'/■Jh '3 


3re 


y( IT) = y(t=T 


( 6 . 21 ) 


•xmilarly leirom equations (2.50) and (2,51) 


c; / _ 




.Xi+^) + 


^2CT^) 


•) + K,(i ,(T ^ ^- 




Using equations (6,11), (6,3) and (6.4) in above, wc obtain 


X ( T ) 

(Iq) ~ o>^(>C.^ ( T^) — ) - ayClj^) 


(6.22) 


where 


y(Tc) 


d_y(t) j 
dt t=T, 


(6,23) 


From equations (6.2) and (6,13) 


X,( T ) X,(o) 

(y^ir^) -~k— = iMXjCo) --fj-) 


•^ro Xo(o) 

^ (Xi(o) +-Y3— > 


Using oquatioms (6.11) and (6.18) this gives 

(X (X )- i-- (y(o)+ay(o)) - y(o) (6.24) 

2^ C s|3 (2u)„) 


Substituting £^ncom equation (6,24) into equation (6,22) we 


o btain 


Y(t ) = i (1+ y(o) - w^(i+ \) y(o) 

C Z r* (jl) 


(6.25) 
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The equations (6,20) and (6.25) gives the desired boundar 
value conditions required to solve equation (6.10). These 
relate tho end point value and slope of variable y(t) to the 
initial value and slope of this variable. 

6.2.3 Derivation of tho first order equations relating y and 
(dy/dt) 

Let a nov; variable z be defined as 


py = dy/dt = z 
Then this implies that 


P^y = (dy/dt) = dz/dt =(dz/dy) (dy/dt) 

That is. 


P 


2 


y = 


Z 


dz 

dy 


(6.27) 


(6.28) 


Substituting equations (6.27) and (6,28) into equation (6.9) 
vve got 


dz 

dy 


+ 2az 


4 

3 


K,m T (t) , ^ ^ 


(6.29) 


Thus the second order equation (6,9) can bo reduced to two 
first order equations given by equations (6.27) and (6.28). Tho 
boundary value conditions of equations (6.20) and (6.25) become 


y(T ) = ^ (1-r '^^) y(o) + z(o) 

r “ r 

2 

^^^c^ =5(1+ ^^-) z(o) - 

^r 


(6.30) 

(6.31) 
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6.2.4 General properties of the solution of the derived 
aquations 

There are certain properties of the solution y(t) vi/hich 
can be directly noted from the derived differential equations 
(6.10) and (6.29) 

(i) It can be seen from equation (6.10) that if y(t) is a 
solution for torque T^(t), then (F)^'^^ y(tj is also a 
solution for the torque value of F.T^(t). Therefore, from 
equation (6.12) the current should be increased by a factor 
of (F)^^^ to increase the torque by a factor F. 

(ii) If z = f(y) is a solution of equation (6.29) then from 
equation (6.29) it can bo seen that 

dz j 

dy lYo^Zo ryo»“"o 

This property eliminate tho consideration of two quadrants 
of the phase piano plot [l3] because of the symmetry. 

(iii) It is clear from equation (6,10) that if y(t) is a 
solution of this equation then -y(t) is also a solution. 

is the desired torque waveform and has the frequency of six 
times tho invertor frequency (Soc.4.2). Thus this property 
imply that if y(t) is a solution in interval I then -y(t) is 
automatically tho solution in interval IV. 
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(iv) Let us study the behaviour of equation (6,10) at the 
point when the derivative of y(t) is zero, i.e,, py = 0. At 
such a point the equation (6,10) reduces to 


( 2 ^ 2 , 2 . 

(p +a -i-w^)y 


4 lail), 

3 vTtT 


This iraplies that 


2 4 ^l“r ( 2, 2s 

p y " S 'ir™'* “tCtT “ ^"^r^ 


y(t) 


(6.32) 


During interval I, i.(t) is always positive. 

Thus from equation (6.12), y(t) is negative during this interval. 
Thus equation (6.32) implies that d^y/dt^ is always positive 
when dy/dt is zero. Thus, for assumed positive torque function 
T (t), if there is a possible solution of stator current, then 

2 vm 

the form of during interval I is one of the three types of 

those shown in Fig. 5.1. 

6.3 PHASE PLANE ANALYSIS FOR THE CASE OF CONSTA-NT TORQUE 

In this section the procedure to obtain the solution for 
tho stator current which produces a constant torque has been 
given. 

For the case of constant torque, T the equation (6.29) 

qo 

becomes 

dz ^ _ 1 [£ + p2 j _2a 

dy z '-y 


(6.33) 
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where 


a 



K,w T 
i r gc 


(6.34a) 


= (a^-Ho^) (6.34b) 

rho equation (6.63) gives the relationship botwoc-n y and its 
derivative z (= dy/dt). The phase plane plot [l3] can, 
therefore, bo obtained from equation (6.33). 

The slope of the trajectories dz/dy = S say has foilovd.ng 
properties. 

(i) It tends to infinity as z tends to zero or y tends to 

zero. V'/o note that from equation (6.34a) that a > o as we are 

dealing with the case of torque production due to rotor speeds 

loss then synchronous speed and hence T is positive. 

^ c 

(ii) If z‘is very largo ^ cn -2a and is a constant independent 
of y. This is true therefore in all quadrants of y-z plane. 

(iii) Zero slope isoclin^ i.e. M = o 


i = 0 if i = - fe (f f 

This defines a curve in the y-z plane. 

The locus of points in the phase plane wherein the trajectories 
have the same slope is referred to as isocline 
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This curve- has a maxima or minima when 


1=, £L. -j. 0 ^) = 0 

2a ^ y2 ^ 


This gives 


Y = Yo = i 6 


( 6.36) 


:orresponding to this y of equation (6,36) 


z = z. 




za y 


V = Z = + R 

o “■ a P 


(6.37) 


(iv) For the isoclin with slope S t i«G 


dz/dy = S 


From equation (6,33), the equation of the curve in y~z plane for 


the isocline with slope is 


1 / ^ 1 N 

z = - (r + ? y) 

(S+2a) ^ 


(6.38) 


^’his curve has maxima or minima at 


^ « 2 ^ 


(3+2a) y" 


This gives 


(6.39) 
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This Yjjj of equation (6.39) is same as of equation (6,36). 
As equation (6.39) is independent of S, it implies that all the 
isocline curves have a minima at the same value of y given by 
equation (6,39). From the equation (6.38) it is also clear 
that the isocline curves for different values of S is similar, 
except for a different constant multiple for z. 

(v) It can be seen from equation (6,33) that the reversing of 
the signs of z and y, we get the same equation. Thus, the 
phase plane plot need to bo plotted only for 2 quadrants. 

(vi) From equation (6.33) 



For the zero slope isocline, i.e, for dz/dy = 0, at the extreme 
point (Yjjj^;^) from equation (6.36) y^^ = + , Thus at this 

point from equation (6.71) 

= 0 (6.4l) 

dy^ 


Thus the point (YqjZ^) ^ point of inflexion because at this 


point? 


dz_ 

dy 


and 


d^z 

— ^ are zero, 
dy 


Using the above information and S as a parameter it is possible 
to draviT a sot of isoclines throughout the phase plane. These 
are then used to determine the trajectories graphically. The 
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method is shovm in Fig* 6.2, The lines from point P, have 
slopes Sj^ and S 2 which aro indicated by the short line segments 
crossing each isocline. They aro extended until they meet the 
isocline indicating slope ^^ 2 * The point P 2 is midway betvjeon 
the intersections. The process is then continued. The tra- 
jectory is plotted as a smooth curve through the points 
which are thus determined. 


The trajectory satisfying the boundary value relations 

(6.30) and (6.31) gives us the desired stator current waveform. 
To find the stator current satisfying equations (6,30) and 

(6.31) following procedure is adopted. 

Lot some arbitary value y(o) bo assumed for solution. 

Then from equations (6.30) and (6.31), the equation for end 
point is given by 

z(Tc) = ^ (1 + ^^)(y(Tc) “ ^§^)y(o)) + 

T X ^ 

2 

- w^(l+ -^)y(o) 


X ♦ G # jt 

z(T^) = Y(T^) + ivl 2 y(o) 


(6.42) 


v/horo 


M 


1 


M2 



■¥ 



) 



W ( 1+ ] 

0)1 


(6.43) 


(6.44) 
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Since y{o) is assumed constant, the equation (6,42) is a equa- 
tion of straight line. The intersection of this straight line 
with the trajectory passing through the point y(o) and z(o) 
gives one value of y(T^). This is compared v/ith the y(T^) 
obtained from equation (6,30). In case these are same, the 
assumed z(o) for y(o) choice gives the solution else a new 
trajectory on the line y(o) = constant is followed. 

The phase plane plot with the trajectories for the case 
of Tq^ = 4 Nt m and = 50 rpm has been plotted in Fig. 6,3, 
Tnis has boon plotted only for two quadrants because quadrants 
II and III are similar to quadrants IV and I respectively 
because of property (ii) (See. 6.2.4), 

To obtain the solution for different rotor speed, a 
different phase plane plot has to plotted. It can be shovm 
that with the change of rotor speeds, wo can change the axis 
in such a manner that the isocline plots remain invariant. 
This can bo seen as follows. 

Equation (6,38) gives the equation of isocline v/ith 
slope S. This can be written as 

“ rs+2i) ^y ^ (6.45) 



(6.46) 
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If tho rotor spood is changed, tho phase plane plot is now 

assuined for different torque T given by 

qc now 


u 


■ qc new 


'05 


r ^ow 


(~ 


o 

a“+w^ I 


r 'now^ 




a + 0 ) 


r ' old 


‘qc ' old 


(6.47) 


with such a choice of the right hand side of equation 

remains the same even for different rotor speed but with 
new value of constant torque. So, if the now z-axis scale is 
assumed as (z/a), the isocline plot will remain the same. It 
should be noted that since the scale of z axis is changed, the 
value of S for each isocline will change but the plot of iso- 
clinnos need not bo changed. 


The phase plane solution gives the values of y(t) as a 
function of y. The time interval can be evaluated from this 
information. The time interval, t^ between y(o) and y^ is giver 
by 


^1 

ti = / (l/z)dy (6.48) 

y(o) 


Thus from this information of y(t) as a function of time, the 
stator current ’waveform can bo evaluated. 


Since y(T ) point corresponds to the end of interval, 

thus from equation (6.48) 
y(t^) 

Tc = , (l/z)dy 

^ y(o) 


(6.49) 
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A interval corresponds to a 60® interval. Thus inverter 
frequency, w, is given as 


c 


(6.50) 


Thus , 


•n: 


CO 


““yTfp”— 

3[ / (l/z)dy] 

y(o) 


(6,51) 


6.4 GEI\’ER.4L NUrAERICAL SOLUTION FOR THE STATOR CURRENT 

In this section, the equation (5.10) has been solved to 
obtain the stator current waveform v/hich produces the 
torque waveform T (t). For the case of stationary rotor, an 
analytical solution can be obtained for the stator current 
waveform which will produce the torque T^(t), For the case 
of rotating rotor, the stator current has been obtained by the 
numerical solution of equations (6,29) and (6.27). These equa- 
tions are the first order equations corresponding to equation 
( 6 . 10 ) . 

6,4.1 Case of stationary rotor 

For the case of stationary rotor, = O) , the equation 
(6,10) reduces as 

(p^+2ap+a^) y(t) =0 


(6.52) 
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For the solution of equation (6.52) these are tvi?o boundary 
value conditions. These- have been obtained by substituting 
= C in equations (6.18) and (6.22) ivhich gives 

y(o) = -ay(o) (6.53) 

y(T^)= -ay(T^) (6.54) 

The S'jlution of equation (6.52) is 

y(t) = Y^e^^^ + Y^t (6.55) 

v/hexe Y^^ and Y 2 are constants. These are evaluated using bound- 
ary value conditions. Substituting from equations (6.55) into 
(6,53) vjQ obtain 

-aY^ + ^2 = -aY;^ 

i.e., ^2 (6.56) 

The equation (6,54) gives the condition oiniulati<=>R. as (6.56;. 
Thus, from equation (6.12) 

Mt) = - (-irV V*’ 

" cl ^ 

The equation (6.57) gives the form of stator current which 
produces the torque waveform T^(t), 

In case Tq(t) is assumed to bo constant, it can be seen 
that equation (6.57) reduces to the form of stator phase current 
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obtained in See, 5 o i i 

• 0 . 1 , 1 , 

torquG . 


which is shown to produce a constant 


6.4.2 Case of rotor 

The solution of tho equation (6.29) satisfying the 
boundary value coneJi-j- given by equations (6,30) and (6.31) 
Cano bo obtained through numerical analysis. It should be 
noood hole thau hero T (t) is assumed to bo some- arbitairy 

M 

dosirod function. 


Lot the variable y y gf some arbitary time t = t, 

o o 

with r = Zg, Then from equation (6,29), 


dz 

'dy 


y^Yo 



(a2H-4)yJ 


2a 


(6.58) 


To tho first dogroo of approximation, the values of time 
and variable z at y (y^^/^y)^ will be given as 




y=yo 


dy 


(6.59) 


''|y= (Yo+'^y) " '^0+ 

where dy is a small change in y from the value of y^. 

Tho stops for solving equation (6.29) satisfying oqua- 
tions (6.30) and (6.31)^ g assumed constant rotor speed 
are given as follows. 
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» 

(i) Arbitary values of y(o) and r(o) arc assumed for solving 
equation (6.39) 

(ii) y(T ), and z(T ), the end point values of the solution 
of equation (6.29) are calculated from equations (6.30) and 
(6,31) respectively, for the assumed values of y(o) and z(o). 

(iii) The equation (6.29) is solved from y(o) to y(T^) in stops 
of i} y, using equations (6.58) and (6.60) and the value of 
z(T^) is obtained. 

(iv) Values of z(T ) obtained in steps, (ii) and (iii) above 
are compared and the difference is calculated, 

(v) Stops (ii) to (v) are repeated with a modified value of 
keeping Yq as in stop (i), till the difference in step 

(iv) is reduced below a limit. 

(vi) Stop (iii) compute the time t, when the value of y has 
has become y(T^), This time t corresponds to the value 
From this, the inverter frequency, can be computed, as 
from equation (6.50) 

w = 7 c/3T^ (6.61) 

The step (iii) gives the solution of equation (6,29) 

'When the assumed initial values of y(o) and z(o), give the 
difference of stop (iv) less then a limit. From this, the 
stator current v/aveforr*! during this interval can be obtained 
using equation (6,12). 
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Having obtained the stator current v/aveform for an 
interval of inverter, v/hich is inverter I in our case, the 
stator current for other intervals can be obtained using the 
properties given in Sec. 5.2. 


This stator current will produce the assumed torque wave- 
form T ( t) . 

q 

6.5 NUi'ffiRICAL SOLUTION FOR THE CASE OF CONST/UNT TORQUE 


Equation (6,29) has been solved for tho case of constap-t 
torquo. This equation for T^lt) = T^^ becomes 


^ dz , o „ 
z ^ + 2az 


4 /^l^r^qc 

3 hi 


) 


/ 2 2 \ 
(a + 0 ) )Y 


(6.62) 


If a now variables are defined as 

■3{- 1 -K- 

y , z - 

a 



(6,63) 


Then from equation (6.12), equation (6.62) becomes 


dz , o * 
z + 2az = 

dy 


4 1, 

Q \ T /4s 

qc y 


/ 2^ 2. 

-* (a +w ) y 


(6,64) 


* -Sfr 

Tho boundary conditions in terms of y and z are similar to 

4C" ^ 

equations (6,30) and (6,31) with y and z replaced by y and z 
respectively. 


This equation is solved using stops given in Section 6.4 
through a computer program, whoso listing is given in the 
Appendix B at various values of rotor speeds. 
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Figs. 6.4 to Fig. 6,9 give the plots of initial points of 
solution and corresponding inverter frequency for various value 
of the rotor speeds. From these plots it can seen that 

(i) There is a lower limit of frequency in the plots. This is 
so because these plots have been drawn for a particular rotor 
speed and the constant torque is assumed to be positive. 

Thus the synchronous speed corresponding to the inverter fre- 
quency has always to be more then the rotor speed, 

(ii) The plot of inverter frequency with respect to the initial 

jf 

value of the solution of y(t) is parabola. This parabola opens 
up as the rotor speed decreases. In case of parabolic form, 
there is a upper limit of frequency. 

(iii) There a lower limit on the value of y(o). This is so 
because protection of a constant value of torque has been 
assumed. The current is inversely proportional to y(t). The 
initial, slope of y(t) is negative and so the slope for the 
solution is always negative (Sec, 6.4.2 point (ii)). Thus the 
lower limit on y(o) corresponds such a initial current vhich 
M'illl produce torque value greater then assumed torque value. 

Thus there is no solution during this region. 

(iv) For a particular value of inverter frequency there is only 
one solution for stator current which will give us the assumed 
constant torque. 
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(v) Points of interest of lov.; slip one on the lower branch of 
the parabola. On this as the initial value y(o) value in the 
solution increases, z(o) value becomes more negative and the 
corresponding value of the invertor frequency decreases. 

(vi) From the Fig. 6,4 to Fig. 6.9 it is possible to plot the 
initial values (y(o), 2 ( 0 )) for the solution for the various 
jof rotor speeds. This has been drawn for the inverter frequenc 
of 20 Hz in Fig. 6.10, for the constant torque at all rotor 
speeds. This plot can be used to obtain the stator current 
waveform to obtain constant torque at various rotor speeds for 
the inverter operating at certain constant frequency. 

Fig. 6.11 compares the solution of nonlinear equation 
at 500 rpm to produce a constant torque and the optimum exponen- 
tial waveform for this torque and rotor speed [sec. 5.3], It 
is seen that there is a close correspondence. This implies that 
exponential waveform corresponds closely to the stator current 
required to produce constant torque. 
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sp<?ed = tOOrpm 
torque = 6 N-m 
y’= I l<a 




Fig 65; InitialT Rjint Of The Solution 
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CH/J^TIZR 7 

COR'CLUSICNS 

The Dnalytical expreessions for the reforrod rotor current; 
of the induction r.otor fed by ideal square wave current have 
been obtained. The rotor current has been computed by solving 
motor performanco ocur'tions in dq frame. They are also comput- 
in';, by calculating the transfer function of the motor. It has 
been shown that both these methods give the same results. 

The analytical expressions for the induction motor being 
fod by stator current with cosinusoidal rise and fall during 
aoiomutation have also been obtained. A more complex time 
domain expressions are given. It has been shown that as 
w “p this solution approach the solution obtained for 
sc;uare wave current. 

Tho nature of tha electromagnetic torque produced by 
the motor has boon studied. The general expression for the 
torcuo in a reference frame revolving at arbitary speed has 
been given. It has boon shewn that the torque is independent 
of the speed of the frame. It has boen shown that the torquo 
harmonic have frequencies of the type Rpm, where w is the 
frequency of stator current and p is a integer. For the case 
of motor being fod by an invertor, the torque frequencies are 
of the type 6pio. 
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The methods to compute electromagnetic torque have been 
given. Tliis methods are applicable to both ideal and nonideal 
current source inverter. This torque has harmonics (Fig, 4,5). 
It has been sho'.vn that for computation of a torque harmonic of 
the order ’6p' , it is sufficient to consider the current har- 
monics of the order (6p-l) and (6p+l), (Table 4.2). Detailed 
calculation show that interaction of fundamental in rotor and 
those two harmonics in stator contribute the major percentage oj 
the torque. 

The control of the torque hcarmonics by modulation of the dc 
input current has been studied in Chapters 5 and 6. It is showr 
that for the balanced system the modulating frequency should be 
six times the inverter frequency. The cases with exponential 
and cosinusoidal modulations have been studied to obtain con- 
stant torque. The study with exponential modulation shows : 

(i) It is possible to obtain a constant torque for stationary 
rotor with exponential modulation (Table 5.1). 

(ii) It is not possible to obtain constant torque for rotating 
rotor with exponential modulation. 

(iii) At speeds less then synchronous spood, tho exponential 
modulation always improves tho performance. At speeds greater 
then synchronous speed, there is deterioration in the perfor- 
mance with exponential modulation (Figs. 5.2 and 5.3). 
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(iv) At a given rotor speed the parameters of the exponentia 
modulation have a optimum value. With these values the sixth 
harmonic torque can be reduced to nearly z-zero (Fig, 5.4). 

The optimum choice varies with rctor speed (Fig. 5.5). 

"■i'ith cosinusoidal ri2)dulation,also, it is not possible to ge1 

constant torque. Tables 5.2 to 5.6 show that 

(i) There is a distinct phase angle for every choice of 

amplitude and frequency of the cosinusoidal modulating wave- 
form Vi/hich gives maximum performance index 

(ii) = 6w, it is possible to reduce sixth harmonic torque 
nearly equal to zero (Table 5.2), 

(iii) There is not much variation in the average torque with 

With = 3m there is a significant 
variation in the average toirque. 

(iv) For = 3ci> or the4e is an optimum amplitude for 
maximum performance index. For case of = I2w or 18w, 
this optimum nature is not observed. 

(v) There is not much improvement in performance index with 

^m ~ compared to = 6w or 3w, 

(vi) It is possible to reduce 12th harmonic of the torque by 
current with oj„ = 6w (Table 5.3),, 

lu 

(vii) It can be said that in case the cosinusoidal modulation 
is being used, the modulating waveform frequency = 6w gives 




best control. 
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The inverse problem cf computing the profile of the 
inverter output current t.: obteino the desired torque waveform 
has also been studied. The stater current for a desired torque 
waveform is given by the solution of a nonlinear second order 
equation. The expression for current to obtain a desired 
torque 'waveform at stationary rotor has boon obtained. If tho 
desired torque v>favofcrm is alv/ays positive, then the possible 
shapes of the stator current to obtain this torque is given 
by Fig. 6.1. The solution of the second order equation need 
two boundary value conditions. These have been obtained and 
these- turn out to be independent of torque. This second order 
equation is reduced to two first order equations. These can 
bo solved by phase plane analysis or numerical integration. It 
can be- said that given the rotor speed and constant torque value 
it is possible to obtain a current waveshape which will give 
this performance. There is a range c-f invertor frequency in 
Vi/hich this solution is possible (Fig. 6,4 to Fig. 6.9). If 
the invertor frequency is defined in this range then there is 
a unique current waveform. Outside this range, there is no 


s:.'luticn. Tho solution of tho nonlinear equation gives current 
which is close to the exponential waveform (Fig. 6.11). 
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Details of the slip ring induction 


Primary 
4C0 Volts 
5,0 Amporos 


3.0 H.P. 
50 CyclGS 
1400 rpiB 


The parameters of the motor are 

Stator leakage inductance ^11 

Rotor leakage inductance 

roforred to the stator ^ 

Mutual inductance botv/een 

stator and rotor phases 

referred to stator M 

Stator resistance 

Rotor resistance referred to 
stator 


motor used 

Secondary 
145 Volts 
10 Amperes 


= 0.2453 H 
= 0.2453 H 

= 0.2364 H 
= 1.6 Ohms 

= 3.31 Ohms 
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